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SUMMARY 


During  a  previous  program,  Cornell  Aeronautical  Laboratory,  Inc.  , 
developed  a  method  of  computing  rotor-blade  loads  and  motions  of  a  single¬ 
rotor  helicopter  in  steady  forward  flight  by  assuming  that  the  blade  pitch- 
control  settings  (collective,  longitudinal  cyclic,  and  lateral  cyclic)  and 
the  rotor-shaft  tilt  angle  are  known  and  that  the  blade  motions  are 
restricted  to  flapping  and  flapwise  bending.  The  present  effort  was  under¬ 
taken  to  extend  the  previously  developed  method  by  (1)  including  the  blade 
inplane  and  torsional  motions  and  (2)  treating  the  four  trim  constants 
(namely,  the  blade  pitch-control  settings  and  the  rotor-shaft  tilt  angle)  as 
unknowns. 

The  trim  constants  are  different  for  different  flight  conditions  and 
are  determined  through  the  use  of  four  appropriate,  average  equilibrium 
conditions  of  the  helicopter.  These  equilibrium  conditions  are  called  the 
trim  equations  and  are  derived  in  this  report,  taking  into  account  the 
inertial  forces  of  the  blades  due  to  elastic  deformations. 

Lagrange's  equations  for  the  blade  motions  are  given.  The 
generalized  coordinates  employed  in  these  equations  to  represent  blade 
bending  are  those  which  have  deflection  components  in  two  mutually 
perpendicular  directions.  Orthogonality  relations  of  vibration  modes  of 
twisted,  rotating  blades  are  derived  and  are  used  in  simplifying  the 
equations  of  blade  motions. 

A  successive  approximation  procedure  was  developed  which,  upon 
incorporation  with  the  previously  developed  iterative  procedure,  yields 
the  aerodynamic  loads,  the  blade  responses,  and  the  required  trim 
constants. 

Computed  results  are  obtained  for  the  UH-1A  rotor  at  advance 
ratios  of  0.  26  and  0.  08  and  for  the  H-34  rotor  at  advance  ratios  of  0.  29 
and  0.  18.  Comparisons  of  these  results  with  available  measured  results 
are  presented. 
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INTRODUCTION 


The  U.  S.  Army  Aviation  Materiel  Laboratories  (USAAVLABS) 
has  been  directing  a  unified  experimental  and  theoretical  effort 
which  has  for  one  of  its  ultimate  objectives  the  development  of  a  means 
of  predicting  rotor-blade  air  loads,  motions,  and  stresses  in 
helicopters.  That  is,  it  is  desired  to  develop  the  capability  of  predicting 
these  quantities  when  all  the  rotor  hub  and  blade  physical  parameters, 
the  rotor  operating  flight  condition,  the  helicopter  weight  and  center 
of  gravity  location,  and  the  fuselage  aerodynamic  characteristics  are 
specified.  As  a  participant  in  this  Army  effort,  Cornell  Aeronautical 
Laboratory,  Inc.  (CAL), has  been  working  on  the  development  of  a 
method  of  prediction  for  a  single-rotor  helicopter  in  steady  forward 
flight  (Reference  1).  This  earlier  development,  however,  used  the 
assumptions  that  the  blade  motions  are  restricted  to  flapping  and 
flapwise  bending  and  that,  for  each  flight  condition,  the  blade  pitch- 
control  settings  and  the  rotor-shaft  tilt  angle  are  specified  by  four 
given  constants.  In  reality,  these  four  trim  constants  are  unknowns 
to  be  determined  by  the  helicopter  weight  and  center  of  gravity 
location  and  the  fuselage  aerodynamic  characteristics,  as  explained 
later  under  "TRIM  EQUATIONS".  Furthermore,  the  blade  inplane 
(or  chordwise)  and  torsional  motions  could  be  important  (directly  and 
through  their  coupling  with  each  other  and  the  flapping  motions)  with 
respect  to  the  blade  stresses  and  aerodynamic  loads.  Consequently, 
the  present  program  was  undertaken  to  improve  the  previous  develop¬ 
ment  by  removing  the  above-mentioned  assumptions.  The  two 
diagrams  in  Figure  1  show  the  differences  between  the  previous  and 
present  programs. 

The  same  method  and  computer  program  developed  in 
Reference  1  to  predict  the  blade  aerodynamic  load  distributions,  when 
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the  blade  response  and  the  trim  constants  are  given,  are  uaed  here. 
This  means  automatically  that  the  same  wake  and  blade  representa¬ 
tions  (both  described  fully  in  Reference  1)  arc  used  in  both  programs. 

PREVIOUS  PROGRAM  (REF.  1) 


Y  S I  CAL  PROPERTIES: 

HADE  PROFILE  DRAG 
COEFFICIENT 


HELICOPTER  WEIGHT 
AND  C.G.  LOCATION; 

FUSELAGE  LIFT.  DRAG 
ANO  PITCHING  MOMENT 
COEFFICIENTS; 

01  STANCE  BETWEEN  MAIN 
AND  TAIL  ROTOR-SHAFTS 


FLIGHT  CONDITIONS 


BOUND  CIRCULATION 
DISTRIBUTION 


BLADE  LIFT.  DRAG,  AND 
PITCHING  MOMENT 
DISTRIBUTIONS 


GENERALIZED 

AERODYNAMIC 

FORCES 


PITCH-CONTROL  SETTINGS; 
ROTOR  SHAFT  TILT 


AVERAGE  ROTOR  THRUST. 
LONGITUDINAL  FORCE, 

SIDE  FORCE,  TORQUE,  AND 
PITCHING  MOMENT 


BLADE  LIFT.  DRAG,  AND 
PITCHING  MOMENT 
DISTRIBUTIONS; 

BLADE  RESPONSE  AND 
BENDING  MOMENTS 


Figure  I.  COMPARISON  OF  PREVIOUS  AND  PRESENT  PROGRAMS. 

Except  for  minor  modifications  that  are  necessary  for  admitting 
blade  profile  drag  coefficient  as  additional  input  and  for  having 
blade  drag  (profile  plus  induced)  distribution  as  additional  output. 
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TRIM  EQUATIONS 


Consider  a  single-rotor  helicopter  equipped  with  a  tail  rotor 
(which  is  assumed  to  produce  a  side  force  and  yawing  moment).  When 
the  helicopter  is  in  steady  level  flight,  the  flight  condition  specifies  its 
forward  speed,  Vf  ,  and  the  rotor -shaft  rotational  speed,  /l  .  For 
each  given  flight  condition,  the  average  forces  and  moments  (from  both 
the  aerodynamic  and  the  inertia  sources)  acting  on  the  rotor  must  be  in 
equilibrium  with  the  helicopter  weight,  the  aerodynamic  forces  and 
moments  acting  on  the  fuselage,  and  the  side  force  produced  by  the 
tail  rotor.  The  equations  that  express  the  requirements  to  secure  the 
equilibrium  are  called  the  trim  equations.  In  a  forward  flight,  the 
fuselage  ride  force  and  yawing  moment,  as  well  as  the  net  value  between 
the  rotor  average  rolling  moment  and  the  fuselage  rolling  moment,  are 
all  small.  Consequently,  their  effects  on  the  trim  equations  will  be 
ignored.  Then,  besides  the  forces  and  moments  in  the  longitudinal 
plane,  we  need  only  to  consider  the  presence  of  a  tail  force  equal 
and  opposite  to  the  rotor  average  side  force  with  its  moment  about 
the  rotor  shaft  equal  and  opposite  to  the  rotor  average  torque.  Since 
the  consideration  of  the  forces  and  moments  in  the  longitudinal  plane 
leads  to  three  trim  equations,  there  will  be  a  total  of  four  trim 
equations.  With  the  pitch-control  system  contributing  a  blade  pitch 
angle  &  according  to 

e  *  e0  +  o,c  cos  +  e,s  sen  y  ,  (i) 


*  Their  effects  can  actually  be  cancelled  by  secondary  adjustments  to 
the  tail  rotor  force  and  location  (perhaps  combined  with  a  lateral 
tilt  of  rotor  shaft).  In  the  derivation,  the  tail  rotor  location  is 
considered  to  be  fixed. 
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the  constants  &0,  *)lc  ,  &l$  together  with  the  rotor-shaft  tilt  angle,  CtSf 
in  the  longitudinal  plane  may  be  called  the  four  trim  constants.  Their 
values  depend  on  flight  conditions  and  must  be  so  determined  that  they 
make  the  trim  equations  satisfied  for  the  given  flight  condition. 

Referring  to  Figure  2  and  assuming  that  &s  is  sufficiently 
small  that  cos  ots  •  /  and  sen  ots  »  acs,  the  trim  equations  may  be 
written  as 


”0 

2  7T 

rZn 

/  X  d<J> 
Jo 

=  yVas  -  xF 

(2) 

2  7T  ^ 

r2rr 

/  Zd<J>  •  W  -  Zf 
'  0 

(3) 

”5 

2rr  - 

rzjr 

j  M  dtp  = 

o 

(Wa5  -  XP)D ,  -(W-Zf)Dx-Mf 

(4) 

.2ir 

**  1  Yei*  •' 0  • 

n» 

2tt  k. 

/  Qdf  + 

’n 

(5) 

Note  that  is  the  number  of  rotor  blades;  X,  Y  ,Z,  M,  Q  are  force 

and  moment  components  acting  on  one  rotor  blade  in  the  directions 
indicated  in  Figure  2;  Rt  is  the  distance  (in  the  *  -direction)  between 
the  main  rotor  and  tail  rotor  shafts;  and  ne  tail  rotor  produces  a 
force  in  the  y-dircction  equal  to 


For  a  given  flight  condition,  the  blade  loads  (X,  Y  ,  Z  ,  M  and  Q  ) 
depend  on  Ocs  and  the  blade  motions,  which  may  include  flapping 


(  hf  ),  lead-lag  (  ).  pitching  (  6),  flapwise  and  inplane  bending 

deflections  (  h4  ,  j  >/  without  lead-lag  and  j>2  with  lead-lag),  and 
torsional  deflections  (&y). 
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I  C.G.  OF  HELICOPTER 


TOTAL  WEIGHT 
OF  HELICOPTER 


Figure  2.  FORCES  ACTING  ON  HELICOPTER. 
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The  plane  that  is  normal  to  the  rotor  shaft  and  contains  the 
flapping  axes  will  be  chosen  as  the  reference  plane.  Each  bending 
mode,  represented  by  h4  ,  will  be  described  by  two  functions  of  the 
radial  distance,  r  ,  from  the  rotor  shaft  of  a  blade  section:  f, 

i 

describes  the  radial  variation  of  the  deflection  of  the  elastic  axis  nor¬ 
mal  to  the  reference  plane,  positive  upward;  fc.  describes  that 
parallel  to  the  reference  plane,  positive  forward.  (Thus,  for  each 
bending  mode,  h  A  is  the  flapwise  bending  and  h  /  is  the 
associated  inplane  bending.  )  Pure  flapping  or  lead-lag  may  be  con¬ 
sidered  as  a  special  bending  mode  with  -fc  =  O  or  -  0  ,  respectively, 
both  with  zero  curvature. 

In  steady  flight,  the  blade  loads  and  motions  are  periodic  so 
that  they  can  be  resolved  into  harmonic  components.  For  example, 
the  lift  per  blade  may  be  expressed  by 

OO 

L  L0  *21  (Lnc  COS  +  Lns  Si*  n  t)  (6) 

n  m  / 

and  the  generalized  coordinate  representing  the  i**  bending  mode  may 
be  expressed  by 

hi  =  h[l)  +2L  (hnc  COS  nf  +  h(‘l  Sin  n  f).  (7) 

n  *  1 

Note  that  the  pitching,  9  ,  has&nc  =  &ns  =  0  for  n>2.  Throughout  this 
report,  the  subscripts  O,  nc  and  ns  will  only  be  used  to  indicate 
harmonic  components  in  exactly  the  same  manner  as  in  (6)  or  (7). 

Let  J.  be  the  distance  of  local  c,  g.  (i.  e.  ,  the  c.  g.  of  a  blade 
section  at  distance  r  from  the  rotor  shaft)  above  tne  reference  plane. 
Then 

J  *  (r-rc)fic  +Xhjfh  -X/(°d  +  Z  ) 
j  J  j 

-  (e,  -  C  +  X, )(60  *  9,c  cos  i‘  *  &,s  sin  y)  (8 ) 
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where  J5C  is  the  built-in  coning,  rc  is  the  coning  offset  radius  (  /5C  -  0 
for  r  <  rc  ),  9B  is  the  blade  built-in  twist  (about  elastic  axis),  is  the 
distance  from  midchord  point  to  the  pitching  axis,  e  is  the  distance 
from  midchord  point  to  elastic  axis,  and  X ,  is  the  distance  from  c.  g. 
to  elastic  axis.  (All  these  distances  are  measured  chordwise  and 

>jc 

positive  forward;  see  Figure  3.)  The  acceleration  of  the  local  c.  g. 
in  the  ^-direction,  ay,  can  now  be  obtained  by  taking  the  second  deriv¬ 
ative  of  J  with  respect  to  time,  yielding 

*h.  - 
j 

+  (e,  -  e  +  X,) Cl  (9,ccos  *  0,s  simp) .  (9) 

Since  Q-y  is  actually  tilted  toward  the  rotor  shaft  by  an  angle  equal  to 
the  flapping  angle,  h,/(  R ->y),  plus  fic  ,  it  has  a  radial  component 
equal  to 


Thus,  the  total  radial  acceleration  of  the  local  c.  g.  is 

ar'  -(fic  +£r)  -  a‘r  •  <10> 

The  tangential  acceleration  (positive  forward)  of  the  local  c.  g.  is 

«c  *  Z  h,  %  *  n2(x,  -  £  h  f  )  (11) 

J  *  J  J 

where 

/,  *  c,  -  e  v-  xf  -eD  -  e,  (1Z) 

and  e0  and  eP  are  the  chordwise  offsets  of  the  lead-lag  hinge  and  the 
pitching  axis  as  shown  in  Figure  3. 


*  All  the  Coriolis  accelerations  will  be  neglected. 
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NOTES:  (I)  THE  REFERENCE  PLANE  IS  PERPENDICULAR  TO  THE  ROTOR  SHAFT 
AND  CONTAINS  THE  FLAPPING  HINGES. 

(2)  THE  BLADE  REF ERENCr  LINE  LIES  IN  TH_  REFERENCE  PLANE, 
INTERSECTS  the  ROT  "  SHAFT,  AND  IS  PERPENDICULAR  TO 
THE  FLAPPING  HINGI 

Figure  3.  BLADE  Ot  FSETS  (POSITIVE  SHOWN). 
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Letting  m  be  the  blade  mass  per  unit  span,  -l  be  the  lift  per 


unit  span,  dp  be  the  profile  drag  per  unit  span,  and  t  be  the  induced 
angle  (so  defined  that  the  total  drag  per  unit  span  is  equal  to  dp  t  ej.  ), 
the  contribution  per  unit  blade  span  to  /  can  now  be  expressed  as 


—  *  (dp+ €  1  *  mac)scn<t-  )  2  *  mcir  J  cos  + 

-  Idp+eJ  +m]T/i.{e  +mn.l(x3  -£hjfc 

L  j  J  j  J  '  > 

■  f(A  ''”*>)  -  mnrjcos  . 

(13) 

Similarly, 

dr  "{dr  +  *2  *  mY^hi  \  +  ma  (*3  *Z  h,  %  )}  co5<P 

J  J  V  j  J  ,J 

*  l(^c  *R~^7)(J'rna'>} ' mfl  r\ 3in  * 

(14) 

1 

II 

M  v 

(15) 

*7  ■  (d,  +  '  "L  v.fe. )  r 

(16) 

and,  letting  I  be  the  blade  mass  moment  of  inertia  (about  the  elastic 

axis)  per  unit  span  and  be  the  blade  aerodynamic  pitching  moment 

(about  the  elastic  axis,  positive  nose  up)  per  unit  span, 

~~  -  -wa-j )  r  -  njflrJ  j  cos 

S  +(‘*r  +  €'t)[(r‘rc)Ac  +  Z-  hj hr  (e< -e)e} 

J 

+  (ja'm  X,  Z)(rfd<c  cos  V  *  n*&f3  sen  <P  ~Z  te  ) 

+  -(e,  -  e  -eD  -eP) 

(17) 

9 


.  2rr 


In  obtaining  /  X  dp  .  ....  /  Mdp  from  Equation ( 1 3), 

•'o 

....  Equation  (17)  by  double  integration,  it  is  essential  to  notice 
the  following  relations: 

,  2rr  ,  2rr 

2rr  ; 


J  dp  *•  2rr  ;  J  cos2 p  -dp  -  J  s<.n*p-dp 

/Zrr  ,  2rr 

cos  np- dp  =  J  sin  np  -dtp  *  O  ,  r>  >  f 

(irr, 

J  (cos  or  stnnp)  •  ( cos  p  or  sm  p)  •  dp  ■  O  ,  n  >2 

(^j)o  m  (&j)o  =  '  °  >  (a'})fC  =  ■  ^  ire 

~n*h'c  ■>  (*j),s  • 

&),C  -  -rfelic  >  (6j),s  =  -nle,(sJ> 


For  example,  taking  the  advantage  of  the  above  relations,  the  con¬ 
tribution  of  the  first  line  of  Equation  (17)  to  J A/^can  be  obtained 
as  follows: 

ru>  -  ma.^)  r  -  mCi  r  $  |  cos  p  •  drj  d  </ 

=  J  [  J  {dte  -  m(a.})  tc  -  mn*$fC]^cos*p-dp^r  Jr 

f* 

/  ^fc  rdr  r  ttB,c 

J  n 


-  V 


where  B,e  is  the  amplitude  of  the  first  harmonic  cosine  component 

/# 

Jrdr  ,  The  actual  procedure  of  computing  3,c  involves 
calculating  B  at  2N  equally  spaced  azimuth  positions.  Denoting 


k-t 

/V 


TT 


and  0*  w 


(k  •  1.  2.  ■  ,  2 N)  ,  &  ,c  is  given  by 


2N 


' !C 


~  3k  cos 


N 


k-t 


See  Appendix  III  of  Reference  1. 


10 


If  the  average  value,  B0  ,  and  the  amplitude  of  the  first  harmonic 

>!< 

sine  component,  &fS  ,  are  also  wanted,  they  can  be  computed  from 

8o’2N^—Bh  and  3>s’~y  aksc"i>k. 

k  *  /  k-r 

Zrf 

The  contribution  of  the  second  line  of  Equation  (17)  to  /  Mdf  can 

Jo 

be  obtained  by  writing  it  as 

(fac)ti  jo  *(^e)o  fta  |  dr 

-  ^ 

-  rr  Jo  (dp  *  el)  „  [(r-re)fic  +y  h[j)f X.  -  (e,  -  e)  dr 

-  7 t (dp+  ej,)0  h{s>Sf,j  -  (ef  -e)&f5^dr  . 

However,  in  the  first  integral,  an  approximate  expression, 

J  •  (r  -rc  )  ht  ~  (*’  -  g  *  X,  )  & , 

will  be  used  instead  of  Equation  (8)  to  make 

( d-c)is  j-O  +  (O'c)o  }f3 

•  •  in ‘ ( Z. h<ts % ){(>"•  rc) Ac  *  hS?-(t,-e+J.,)5, j 

*  n*  (x,  -  £/.'%)  { hi?  -U,-e*7,)  §„  } 

while  in  the  second  and  third  integrals  the  replacements  of 


and 


with 

j  J 


r  -  rF 
R-rF 


2  h  -fh  with 
</  J 


r-  rF 

R-r, 


h 


(t) 

is 


will  be  made.  The  complete  set  of  results  obtained  by  double  inte¬ 
gration  of  Equations  (13)  through  (17)  is  as  follows: 


*See  Appendix  III  of  Reference  1. 
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ij  W'-OH'Zn'ZmC'KK-^-l-o 

Vq  j  i  a  r f 

ipJ^Zdt--  2t, 

1  /■*■ 

wj0  <**  ’  ^ 

i-f  W-  B,c  *ZSX*£h%' \Ac(^  -r,mCj) 

^  J 

h(f)  /  )  h<r) 

*-r£(,vr"% ) -a‘<rt  ”<"'>*,  +t*9»*-'.’. 

+  at£M°f  l  J^r  (*rcj  “  '>  "’CJ  )  -  0fS  T*t  CJ  j 
•  U(n  h(° 

-A  Cn  -^r  C"  *  4  C  ~^fr  c'  *  en  c;" 

-ale„ju  -n'JTe^  t*  -(Mp) 


where 


X2  =  e,  -  e  *■  X, 


m‘i  • 

1  m%  dr 

0 

^2 

r*A  ' 

f  mX,f^  dr 

7*2 

Jr0 

J  mrfc  dr 

ro 

Z1  C 

*1  QJ 

II 

«>? 

XV 

j  mrXi  dr 

H(r-rf)1t 

J  - 
*2*! 

Jo 

r* _ 

j  mX2X3  dr 
r0 

'pi  IS 

X3  *  *z  -ec  -  eP 
[* 

/  m  f*  dr 


dr 


r6 


(19) 

(20) 

(21) 


(22) 
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r* 

“  J  (Joc  m  X?)  dr 

Jro 

L  -  /  UUr 

“ o 


B 


R 

X  rdr 


C '  *  J  (dp+  e  2)(r  rF)  dr 

rr 

c' “  J  (dp*-  el)(r  -  rc)  dr 


rR 

Jgj  -  J  (la~  mX*)  f0jdr 

'Jro 

D  ■  J\dp  +  el) dr 

f* 

C  m  J  ( d-p  +  6l)  r  dr 


C"*f  (dp  +  el)(e0-e)dr 
Jre 

Mpm  Jo  'e  'e0  -er)j]dr 


Note  that  the  lower  limits  of  the  above  integrals  are  based  on:  1^=0 
for  ri  ry  ,  fc.  =  0  for  rsr#(  5  =  0  for  r  <  r0  ,  r>  <  rD  <  re  ,  and 

rc  «  ro  •  For  a  teetering  rotor,  rF  -  r0  =  0.  To  be  strictly  correct, 
the  Za  in  Equation  (18)  or  (19)  must  be  based  on 


the  Lfc  in  Equation  (18)  or  L,s 


when  it  multiplies  J&c  and  on 


’I'^r 


in  Equation  (19)  must  be  based  on 


L 


R 

Udr 


when  it  multiplies 
must  be  given  by 


when  it  multiplies  fic  . 


L  -  f*Jdr 
Jrr 

;  and  m^,  also  in  Equations  (18)  and  (19), 

mhj  *  J  » 

rC  J 

(These  corrections  are  usually  either  not  needed 
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or  are  negligible.  ) 

Define 

Y  =  f^Yd* ,  etc. 
TT  J 0  Tr  'Jo 

Then,  Equations  (2)  through  (5)  become 


X  =  /  (JVccs  -X,) 
n8 

Z  -  ^  (W-Z,) 

M  -  £  \Was-Xp)Dy  ~(W-Z,)Dx-*f,} 
—  —  Q 

r  -  Y  *  e]  s  0  » 


and  Equations  (18)  through  (22)  become 

h">  .  . 


(23) 

(24) 
(25' 

(26) 


A  v*  a  • r**  Ml 


Z  "  2-Lt 


M  -  B,c  pc  C'ri  +  60  C',s  +  6,s  C*  -  ff[r^  (po  is  *  hs  Co  )  ' ^p)is 

*  JTPr  H(r +X-  e<‘  | 

rXlCj  -  pt  (r,c>  -r,ne.) j 

J  ill)  1 

-  2n‘Z_  !>n  [flc&rCj -\mct)-eo**lcl  fp7r(z"i 


T*  'D, 


( f )  Af)  f 

■°«  -irh  L°  -  (fi‘  *A)  l“  * n*  {*£”«,  *« 

-  (fie  f  ipr  )(Z  mk:,  h<t‘  r*, »,  8”’  -  T*,8")}  *R-'  L‘ 

’  J  J 
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The  numerical  values  of  the  right-hand  sides  of  Equations  (27)  through 
(30)  depend  on  the  computed  blade  air  loads  and  motions;  they  will  agree 
with  the  desired  values  given  by  Equations  (23)  through  (26)  only  when  the 
correct  set  of  trim  constants  (as ,  &0l  6iC  t diS)  is  used  in  the  whole  computa¬ 
tion.  The  intricate  dependence  of  blade  loads  on  trim  constants  makes  the 
determination  of  the  trim  constants  for  a  given  flight  condition  difficult. 
Only  methods  involving  iteration  or  successive  approximations  can  be 
used.  One  of  such  methods  is  described  in  the  following  paragraph. 

Let  ( Ocs ,  &Q ,  9fC,  0,s)  be  an  arbitrary  set  of  trim  constants, 
different  from  the  correct  set  (  Ors,60  ,Olc,  &ts)  according  to 

*  <*s  ~  ocs 

&o  *  0o  -  Bo 

B,c  *  0,c  ~  K 

Bft  ~  0fs  ~  Bfs 

When  the  arbitrary  trim  constants  and  their  resultant  computed  blade 
air  loads  and  motions  are  used  in  evaluating  the  right-hand  sides  of 
Equations  (27)  through  (30),  the  resultant  values  will  be  designated 

K,Z,M  and  T 

accordingly.  In  a  similar  manner,  the  following  designations  will  be 
made: 

*a  ,  Z*  ,  ^  .  and  Ta  pertinent  to  (  Ocs  +  As  ,  0O  ,  &fc,  9,s 

Xt,  Zb  .  .  and  fb  pertinent  to  (  as  ,  60  +  A0 , 9ie,  0IS)\ 

Xc  ,  Z. e  .  and  fc  pertinent  to  (  as ,  6 0  ,  9,c  +  Af  ,  9,$  );  and 

;  *4  ,  *  and  Td  pertinent  to  {  as  ,  90  ,  9,e  ,  Bt $  +  A*). 

With  reasonably  chosen  A's  and  (  0cs  ,  0O , 9tc,  9,a),  it  is  reasonable  to 
assume  that  the  following  four  equations  will  determine  a  set  of 
estimated  values  of  as,  90  ,  Gtc  ,  and  §ti  to  yield  a  new  (  cti ,  0O , 6,c  ,  9,s  ) 
set  closer  to  the  correct  set. 

X  +  *  a.  &o  *  *c  @ic  +  &/s  *  X  (2 1 
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*4  S'i  =  * 


L  ^  +  Zb0o  +  Z  '  §te 
&  +  *  *b  60  *  Mc9,c  +  A/d9,s  =  M 

T  +  Ta  as  *-  Tb  90  *-  rc‘  §,c  *■  Td  9/s  =  O  , 


(33) 

(34) 

(35) 


where 


(  )a.  = 

(  >;  - 


()»-() 

As 

C)b  -r ) 


(  )' -  j 
r  AfC 


A 


ts 


(  )  =  X.Z.M, 


or 


T. 


(36) 


The  quantities  X,  Z  ,  and  M  appearing  in  Equations  (32),  (33),  and  (34) 
must  be  computed  from  (23),  (24),  and  (25).  The  aerodynamic  loads  on 
the  fuselage  for  a3  -  ars+  as,  \ocs\  small,  may  be  expressed  by 


Xfi 

- 

X 

ocs 

Zf 

* 

*  ZF 

<xs 

MF 

/V 

■»  M F 

+  Mp 

(37) 

A  / 

where  Xf  ,  XF  ,  etc.  are  constants  depending  on  flight  conditions  and  the 
chosen  reference  shaft  angle,  <Xf  .  Substituting  (37)  into  (23)  through  (25), 

Z 


X  '-^(Wccs  -*r  -X'f2cs) 

nB  v 


M 


*  4"  {(Wa'5  -Xr)Df-(w-Zr)  Dx  -  MF  -  (XF  D}-ZF  +  MF)  &s)  . 


With  0ct-  ocs  +  as  and  Ocs  =  - ( ocs  -  Cls  )  +  Cts  ,  the  above  three  equations 
may  be  written 

X  x  —  i  Was  -  XF  -(W-XF  )  (cis  ~  &s  )  -h  —tW  -Xf  )  ccs 
na  t  *  nB 


(38) 


16 


Z  =  lw -ZF  +-Z'r  (a6  -  S.s)  J  -  ~  Z'F  ocs  (39) 

&m  £  (Was-X^Dj  -(W-ZF)Dz-Mr  -{(W-X,  )Dy  +  Z'FDz  -  AfFj  (as  -  a5) 

+  4  \w~Xr)  +  Z'FDz  -  Mr  j  a5  .  (40) 

By  substituting  these  into  (32)  through  (34)  and  rearranging,  (32) 
through  (35)  become 

<-{xZF-(*-X$as  ■  ffwar^-tW-x^as-ciA-x  (41) 

+Zc&tc  +Z,4&,J  +(za  +  ^~  z'r)as  =  “  jw-Z^  *-z'F  (ac-s  -<x3  -z  (42) 

<  ft  *  xj' 4  /  a/;  ft,  ^  -4.  [i {w-x'F)o}  +zfdz-  m; |J  <*s 

=  -( W-ZF)  \w-Y'f)D}  +ZF  Dz  -MF\(as-as)]-M  (43) 

+  Tc&fc  *■  TJ&rs  +  T^oLs  =  -T.  (44) 

These  are  a  set  of  four  linear  equations,  from  which  the  four  unknowns 
(  as,  60 ,9tc  »  and  0/s)  pertaining  to  a  given  (  ccs  ,  0O  ,  Gtc  ,  9,s  )  set 
can  be  solved.  A  new  {  (X3 , 60  ,  9,c  ,  0/s  )  set  can  then  be  formed  in 
accordance  with  (31)  and  the  computation  can  be  repeated  until  &s  , 

90  ,  6fc,  and  ^become  sufficiently  small. 


I 
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EQUATIONS  OF  BLADE  MOTIONS 


When  those  Coriolis  forces  which  do  not  have  IT  as  a  factor  are 

neglected,  Lagrange's  equation  of  blade  motions  with  respect  to  any  ith 

$ 

generalized  coordinate  (  )  has  the  following  general  form: 


(45) 


j 

where 


%  or  0j  )*(/>,.-  •  •  ,  hs  •  e, ,  •  •  •  • ,  er)  • 

There  will  be  a  total  number  of  (  5  +  /  )  equations  of  the  form  (45). 
According  to  Reference  2,  ^  is  called  the  generalized  mass  coef¬ 

ficient,  g.  the  generalized  Coriolis  force  coefficient,  T„  „  the 

“i 

generalized  centrifugal  force  coefficient,  and  Kg  the  generalized 
restraining  (elastic)  force  coefficient.  (Formulas  for  these  coefficients 
can  be  obtained  by  the  method  given  in  Reference  3.  )  Symbolically, 

indicates  the  presence  of  a  generalized  structural  damping  force 
(associated  with  y  )  equal  to 


h  ^  ? 


OO  . 

J~  Stn  n  ^  +  q  (j)  cos  n 


n  :  t 


ns 


*) 


(46) 


7, 


H  (9^  C°S  ”9  +  9ns 


n  z  i 


(47) 


The  constant  for  each  ^  is  called  the  structural  damping  factor  of 
the  degree  of  freedom  represented  by  ^  .  6y  is  the  generalized  aero¬ 
dynamic  force  in  coordinate  ^  and  can  be  calculated  from  one  of  the 


jJj 

If  the  blade  pitch-control  system  were  rigid,  it  would  contr.bute  a 
blade  [  itcti  angle  6  given  by 

G  -  eo  -h  e,c  cos  ^  ■+  6ti  Stn  <p. 

For  a  flexible  control  system,  the  equation  of  motion  with  respect  to 
pitch  deflection  (which  is  the  difference  between  9  and  the  actual  pitch; 
is  similar  to  that  with  respect  to  a  torsional  mode  9i  having  =  1. 


18 


following  two  formulas: 


^  '  JM^dr  -  JUp^i)fc/r 

~  J  ty  ^6.  d-r 


(48) 

(49) 


Finally,  is  the  generalized  force  in  ^  when  the  blade  is  rotating  in 
vacuum  with  all  h  and  0  held  zero.  Cl  and  C#  depend  on  the  blade 
configuration  (  Pe  ,  and  offsets)  and  pitch  control  (  0  ),  and  can  be 

5|t 

calculated  from  the  following  two  general  formulas: 

Ch  -  Ci  Q0  Jm  r  X  ,f^^dr  -rffir  Jmri^  dr 

±Cl(6,c  cos  f  +  9,ssin  </){J mrX,f^.dr  -  Jm(c,  -e  +Xf)S^drJ 

•rfjmt-e,  -eP  ■*■£,  -e  +X,)fe  dr  +  ffjmr X,&a-fh  dr +- rf JmrX,fc.dr  (50) 

(-  eD  -  eP  -  ejj  fg.dr  +  Cl  (6 tc  cos  ^  /  9tssm  <t)(eD  +eF)J/rtX,  f6  dr 

+- Cl  ficJ/nrX)  dr  - 0? J mrdrj-f^  dr 

-CL  f^TatmX,(-ep-eprc,  -<?)J  Ggf^dr  *51* 

where  kA  is  the  radius  of  gyration  of  the  local  blade  sectional  area 
which  is  effective  in  taking  the  centrifugal  force.  Due  to  the  blade 
built-in  twist,  6g  ,  the  centrifugal  force  is,  in  effect,  accompanied  by 
a  torsional  moment  equal  to 

-n‘(fr  mrdr)  ^4  k\ 

*  Unless  noted  otherwise,  a  prime  indicates  derivative  with  respect 
to  blade  radius.  The  first  and  second  X,  appearing  in  (50)  must  be 
replaced  with  (e,  -  e  +  X,  ),  if  h ;  represents  pure  flapping;  the  last  X, 
appearing  in  (50)  must  be  replaced  with  (- ep  +  Ct  -  c  +  X ,  ),  if 
represents  pure  lead-lag. 


I 
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acting  on  a  thin  blade  section  at  r  .  The  contributions  to  Ch  and  Cq 
due  to 

0  ■  0O  e  Olc  cos  4'  *■  &ts  Sin  y 


may  be  first  written  as 


and 


where 


-  Cl2  Th^8) 

~  ( #  8  -  Cl2  T0  #  &  )  ,  respectively, 

Th  e  j  dr 

*  -J m(e,  -e  *-X1)fh.dr 
«  ~/{T<x  +  **><,( -e0  -eP  +  e,  dr 

T&t  g  r  e  *  (&p  +  J m  X  if  dr . 


It  is  now  easy  to  verify  (50)  and  (51). 

(Cfi  U  =  °  for  ^2. 


It  may  be  noted  that  (  C0  )  - 


If  viscous  damping  exists  in  a  certain  jth  degree  of  freedom 
(such  as  lead-lag)  and  does  virtual  work  equal  to  ^  to  a  certain 
degree  of  freedom,  only  the  term  Dm  /n  needs  to  be  added  to  2  Q  for 
the  appropriate  set  of  (  c  ,  j  ).  The  part  of  2  (S’.  that  really  represents 

fi  fj 

Coriolis  force  is  usually  negligible  for  blade  loads  computation.  For 
a  blade  with  offsets  shown  in  Figure  3  and  with  a  rigid  link  connecting 

the  flapping  and  lead-lag  hinges,  the  following  formulas  for  the  M-, 

* 

r-,  and  K -coefficients  have  been  obtained: 


h  - 

J  m  (^hi  ■  fc  ) 

\  dr  +  Jl 

i  j 

Jr  \  1  J  *  j/ 

rD  * 

9  = 

-  f  mX,-fh  dr  = 

Mgjhi 

*  J 

%  ‘  J 

J  A 

h 

r-  rD 


The  lower  limits  of  the  integrals  are  based  on  the  assumption  that 
r0*rD>rfz  0.  (Note  that  6  -  /<>.  =0  for  rsrg  .  ) 
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h 

J 

-  /  (*h 

JrD  •  J  1 

£  *  fcl  \  )(Jr  mrdr)dr 

-  |  Jh  -  2  Jut  +(r 0-^)Jr  mrdr 

+  rf  (r*  - 

<-,)  m„}  (fhi  fi  )rir° 

Th 

« 

9 

J 

'  f  rnrXt-f'h  -f6.  dr  - 
Jrg  'l 

T»i\ 

’  ~  Me  e  -  mX,(-eD 

‘  J  Jr0 

-ep  re,  -  cQ)  fet  fej  dr 

l 

h 

J 

f'  Si.’?  ~  cos  <9  -  fc  ‘ 

i-  tlc  (fc  '  co  s  O  +  -fh  stn  @)[fc  cos  6  *■  />,  sin  6^ dr 


K 


h  e 
c  j 


Ke<*j 


=  O  =  K&  h 
j  1 

=  J  j^GJ  +  il2 ^  j  mrdrj  j  -f'e  dr 


where  ETB,  EIC  and  GJ  are  the  elastic  properties  of  the  blade  section 
(flatwise  flexural,  chordwise  flexural,  and  conventional  torsional 
rigidities,  respectively);  JH  ,  JHZ  ,  and  are  the  inertial  prop¬ 
erties  of  the  link  connecting  the  flapping  and  lead-lag  hinges  (  JH  is 
the  moment  of  inertia  about  the  flapping  hinge,  J**  is  the  part  of  Jn 
due  to  spread  of  mass  in  the  direction  parallel  to  the  lead-lag  hinge, 
r*  is  the  c.  g.  radius ,  and  is  the  mass);  and  &=  §  +  (?g  (but  the 
approximation  6  ~  0O+  &B  has  to  be  used  in  practice).  It  may  be  noted 
that  the  formulas  are  applicable  to  a  teetering  rotor  blade  without 
lead-lag  hinge  by  setting  rp  -  r0  -  0  and  using  appropriate  values  of  JH 
and  (both  belong  to  the  hub  mass  that  flaps  with  the  blade); 
they  are  also  applicable  to  a  rigid  hub  rotor  blade  by  setting 


rf  -  r  q 


=  0. 


One  part  of  the  air  load  program  will  first  compute  each  (r„ 
according  to  the  appropriate  formula  ((48)  or  (49)]  for  ZA  equally 
spaced  azimuthal  positions  and  then,  by  harmonic  analysis,  express 
as  the  following  finite  series: 

n  a/-  / 

% =  (&?Jo + £  +Z  (\  L sinn^  ’ 

r>  */  n  =  / 


it 


(52) 
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which  contains  Zh  coefficients  and  gives  the  actually  computed  values 
of  Cr„  at  those  2N  equally  spaced  azimuthal  positions.  By  making  A/ 

'l 

reasonably  large  (2  N  -  24  for  the  present  program),  it  is  usual  to 
consider  that  it  is  accurate  enough  to  substitute  (52)  into  (45).  By 


doing  this,  all  harmonic  components 
harmonic  are  automatically  nullified. 


By  equating  the  coefficients  of  cosn<p  of  the  two  sides  of  (45), 


£ f  V n‘^< * "  5.  f)\ C 'If',. %  H  9'. 

J 

s(cfiL  *  (\  Lc  ’ 


By  equating  the  coefficients  of  Sin  rnp  of  the  two  sides  of  (45), 

9%  ' n  " 2 5,  1  +  ^  / Kt. s, ' 0 // Ji  ^ 

■  (ct.L  r(atX.  ■ 


Equations  (53)  and  (54)  may  be  called  the  equations  for  nih  harmonic 
components  with  respect  to  .  Note  that  either  (53)  or  (54) 
involves  both  the  cosine  and  sine  amplitudes  of  fa  in  harmonic. 
There  is  only  one  corresponding  equation  for  the  oth  harmonic  compo¬ 
nents  (or  the  constant  terms,  or  the  average  values): 


tX 

*/ 


Since,  according  to  (50)  through  (52),  the  right-hand  sides  ot  both  (53) 
and  (54)  are  identically  zero  for  n  >  N  , 


q(j)  =  =■  o  for  n>N , 

'  nc  '  ns 


(56) 


as  pointed  out  at  the  end  of  the  last  paragraph.  Thus,  there  are 
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(1  +  2/V)  equations  with  respect  to  each  ^  ,  and  each  ^  assumes  the 
form 


?J  i“’  COS  Sinn*).  (57 , 

As  there  will  be  a  total  number  of  ( v5"  +  /  )  different  a  ,  there  will  be 

(5  +  Y)  equations  of  the  form  (55)  to  determine  the  ( 5  +  /  )  different 

aO)  There  will  also  be  2(  5  +  /  )  equations,  half  of  them  of  the  form 

(53)  and  the  other  half  of  the  form  (54)  .  for  each  n  ,  Is  n  <  N,  to 

determine  the  (  5  +  /  )  different  q(j^  plus  the  (  5  +  f  )  different  . 

•  'nc  ns 

Torsional  vibration  modes  of  twisted  rotating  blade,  cantilevered 

at  r  -  ,  will  be  used  for  -f&.  Then  f6,  satisfies  the  following 

equilibrium  condition: 

j^CJ  +  n  *  (J  /rrdr'j  j 

=  uj%.  J  ia /*  dr  -  n2j^  iA  f0.  cosz(e  +  eB)  dr , 

where  u)e  is  the  natural  frequency  associated  with  fg..  Using  cosZlp+Og)**  f, 
the  above  equilibrium  condition  may  be  written  as 

| dr . 

(58) 


Multiplying  both  sides  of  (58)  by  fg.  and  then  integrating  over  r0 
to  R  , 

Jrf  \GJ 

.(u>irO.%’laf,fydr.  (59) 

(Note  that  the  right-hand  side  of  the  above  equation  has  been  integrated 

0R 

by  parts  once,  using  f9  =  0  at  r  =  r0  and  J  Za  0  at  r  -  R  .) 

Using  the  expressions  for  Ke  m  and  .  as  given  before,  (59)  states 

*  j  •  J 

simply 

^6^  0j  ~  (u)&i  . 
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(60) 


Similarly,  starting  with  (58)  for  the  jth  mode, 


*  9-  ~  ‘  )  Met  0j  • 

From  (60)  and  (61),  it  can  be  concluded  that 


(61) 


» 


M0  q  -  0,  j  /  i 

•  J 

(62) 

0^-0.  J  +  t 

(63) 

and 


-  ( u  )0i  '  Tl2)  M  0.  . 

(64) 

Relations  like  (62),  (63),  and  (64)  are  usually  called  orthogonality 
relations  of  vibration  modes. 

Bending  vibration  modes  of  twisted  rotating  blades  will  be  used 
for  and  f  .  The  existing  method  of  Reference  4  yields  modal  shapes 
each  described  by  deflection  components  in  two  mutally  perpendicular 
directions,  from  which  the  and  f  used  herein  can  be  easily  obtained 
for  each  mode.  The  derivations  of  the  orthogonality  relations  among 
these  modes  are  presented  in  the  appendix.  The  results  are 


Mht  hj  -  0  , 

J  *  <■ 

(65) 

h,  *  0  • 

J  *  l 

(66) 

and 

Khihi~-n^Thlhl  *  cuhi  Mhchi  •  (67) 
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By  the  use  of  Km =  2  =  0  and  (65)  through  (67),  (53)  and 

(54)  may  be  written  as 

-nrtzeh^)h£ 

=  (C»t  )nc  *  (&hL  )nc  (<■  -  '  to  5)  (6S 

and 

%(K\*,-  *h>  -nn'te^)h~ 

+  («>\ -"**)  hj*  -*£Z.(»*Mh.ej  +  7*  e.)  e“> 

-  (Chi)ns+(6’\)r,s  (t  -  /  *0  5;.  (69 

By  the  use  of  =  2  ^  ^  =  0  and  (62)  through  (64),  (53)  and  (54)  become 

-  ("*m0i  h.  +  r&i  hJ  h^J  4-  (u>e.  -  »*n')  m0  e  e“c> 

-  n*£  £,  *«'  *  -a')  M.iti  ^  *«» 

(i.  t  to  t)  <™ 

and 

- - A'j ^  e„'; ’  -  n^T 4  *  Ta. h.  j  h'JJ 
*<H  -'’‘n‘)"eieett-rfZi  Tg.g  e<„‘> 

1  j  ->  ‘  J 

=  (fy)**  +  (G-et)ns  (i-t  toy)  (71) 


where 


^  =  r*t*/  *  Met  e- 


r  R  — 

-  7  m  Xt  (-  e  0  -  ep  *■  e,  ~e)  -f0  f0  dr  . 

J  re  1  J 

The  two  equations  for  enharmonic  can  be  obtained  by  substituting 

n  ~  ^ns  ~  ®  int0  (^®)  and  r  -  &ns  -  0  into  (70),  respectively,  and 
replacing  nc  by  0. 


eP  -^e,  ~e)  -f0  f0  dr  . 

*  J 
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The  blade  weight  {my  pounds  per  unit  span)  may  be  taken  into 
account  by  modifying  (48)  and  (49)  as  follows:  subtract  mg  from  the 
first,/  and  add  mgX^o  .  Since  the  blade  lift  usually  amounts  to  over 
20  times  the  blade  outboard  weight,  it  is  permissible  (until  a  more 
accurate  method  of  predicting  becomes  available)  to  treat  the  predicted 
air  loads  as  if  they  included  the  blade  weight  effects.  In  actual 
computation,  the  blade  weight  can  be  more  conveniently  taken  into 
account  by  adding 


' J 


dr 


and 


-h  J  mg  X ,  clr 


to  the  right-hand  sides  of  (50)  and  (51),  respectively. 


During  forward  flight,  the  cyclic  pitch  makes  the  blade  root 
angle  vary  with  the  blade  azimuthal  position.  The  blade  bending 
modes  (and  their  associated  radial  variations  of  bending  moments) 
are  usually  based  on  some  constant  blade  root  angle,  approximately 
equal  to  the  collective  pitch,  In  other  words,  the  blade  elastic  properties 
with  reference  to  a  coordinate  system  fixed  to  the  rotor  hub  are  usually 
assumed  to  be  constant,  approximately  equal  to  the  average  values. 

This  assumption  is  believed  to  be  satisfactory  for  obtaining  the  blade 
motions.  For  the  flatwise  and  chordwise  bending  moments,  however, 
the  following  compensation  may  be  suggested:  after  the  flatwise  and 
chordwise  bending  moments  based  on  the  assumed  blade  root  angle  are 
computed  (by  summing  the  contributions  from  all  the  bending  modes 
employed),  transform  them  to  the  actual  directions  according  to  the 
difference  between  the  actual  (various)  and  the  assumed  (constant) 
root  angles.  (Then  the  resultant  loading  is  taken  by  the  actual  blade, 
although  the  bending  deformation  and  its  associated  inertia  loading 
are  based  on  a  slightly  modified  blade.  ) 
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COMPUTATIONAL  PROCEDURE 


The  trim  constants  are  obtained  by  a  successive  approximation 
method  described  in  the  last  paragraph  of  the  section  entitled  "TRIM 
EQUATIONS".  The  method  starts  with  an  assumed  or  estimated  set 
of  trim  constants  and  then  adjusts  them  systematically  so  as  to  make 
the  average  force  and  moment  components  acting  on  the  rotor  cancel 
the  corresponding  force  and  moment  components  acting  on  the  rest 
of  the  helicopter,  until  the  required  adjustments  are  within  prescribed 
limits.  A  flow  diagram  indicating  the  main  conceptual  elements  of 
this  procedure  is  presented  in  Figure  4. 

The  computational  procedure  described  in  Reference  1  is  for 
computing  blade  loadings  and  response  when  the  trim  constants  are 
given.  It  is  an  iterative  procedure  which  starts  with  a  fir&t  approxi¬ 
mation  to  the  blade  response  and  repeatedly  computes  the  blade 
aerodynamic  loadings  (that  depend  on  the  blade  response)  and  solves 
the  equations  of  blade  motions  (with  the  generalized  aerodynamic 
forces  based  on  the  blade  response  obtained  in  the  previous  iteration 
cycle),  until  the  blade  response  converges  to  within  prescribed  limits. 
The  same  general  procedure  is  used  here  to  compute  the  blade 
loadings  and  response  for  each  prevailing  set  of  trim  constants 
(corresponding  to  the  flow  diagram  shown  in  Figure  4  marked  with 
solid  lines  only).  In  order  to  admit  the  additional  blade  motions  and 
to  include  the  drag  distribution  in  the  blade  loadings,  the  following  four 

modifications  are  made.  First,  the  quasi-steady  bound  circulation, 

$ 

I,  defined  by  (21)  of  Reference  1  as 

,  (72) 


Note  that  Vf  -  flr+  (  l If  cos  cts)  sen  t  and  is  the  slope  of  the  blade 
sectional  lift  coefficient  versus  angle -of-attack  plot. 


i 
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NOTES:  ONCE  «.i(  Sa  .  &,c  AND  £,5  BECOME  WITHIN  PRESCRIBED  LIMITS. 
THE  NEXT  (ALSO  THE  LAST)  CYCLE  WILL  ONLY  ENCOMPASS  THE 
PORTION  MARKED  WITH  SOLID  LINES.  WITH  YES  COMMANDING  THE 
ORAWING  OF  DESIRED  OUTPUTS. 


Figure  4.  MAJOR  STEPS  IN  THE  COMPUTATIONAL  PROCEDURE. 
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has  now 


r 

Oc.  *•  9a  *■  6>0  *■  9tc  cos  p  *  6,s  sin  p  +  2_  &j  f  (* 


(73) 


r  s  s 

Ct  *■  -  rL9,c  s.m  f  +  nets  cob  +  21  &  f  0  +  n(/3  H  (74) 

j*r  J  v  j*r  J  ’ 

s_  .  / 

h  -  Vf  sin  ocs  +  e,  a  (§,c  sm  p  -  6,s  cos  p )  +■  J  hj  fh  ej  &  ie 

,*f  '  j’f  J  J 

3 

*■  (Vf  COS  Ot.s)\£c  +21  h j  )  ros  p  .  (7  5) 


j  *  / 


Second,  the  profile  drag  coefficient,  Cd  ,  of  blade  section  is  a  new 

aP 

input  that  depends  on  Mach  number  (which  is  proportional  to  Vt  )  and 
ae  ,  the  effective  angle  of  at'ack  defined  by 


1  b  . 

ue  =  «+-jae-€ 


where 


t76) 


(77) 


is  the  induced  angle,  with  Ufj/^  being  the  upwash  at  the  three-quarter 
chord  point  induced  by  the  wake.  Third,  the  drag  per  unit  span  is 
equal  to  (  dp  +  Cj,  ),  with 


dp  =  fVl  bCdp  (78) 

being  the  profile  drag  per  unit  span.  It  was  found  in  Reference  1  that 
the  convergence  of  the  procedure  was  improved  by  subtracting  a  part 
of  the  generalized  aerodynamic  force,  approximately  equal  to  the  part 
that  depends  on  the  response  or  generalized  coordinates,  from  both 
sides  of  each  equation  of  blade  motions.  (The  part  subtracted  from 
the  right  or  left  side  is  based  on  blade  response  obtained  in  the  previous 
or  the  present  cycle,  respectively;  they  approach  each  other  as  the 
response  converges.  )  In  effect,  it  means  that  the  convergence  can  be 
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r  t 


accelerated  by  estimating  the  generalized  aerodynamic  forces  for  the 
present  iteration  cycle  closer  than  thos''  based  on  the  blade  response 
obtained  in  the  previous  cycle.  As  the  fourth  modification,  therefore, 
the  terms 


•  h(j) 
h ,  h  "ns 


j*t  ‘ 


e. 


(j) 


ns 


} 


(79) 


(80) 


(81) 


and 


are  subtracted  from  both  sides  of  (88),  (69),  (70),  and  (71), 
respectively,  resulting  in  using  (83)  through  (86)  as  the  formulas  of  the 
equations  of  blade  motions  for  the  present  computation.  A  superscript 
(-1)  appearing  in  (83)  through  (86)  indicates  a  quantity  based  on  blade 
response  obtained  in  the  previous  cycle.  The  formulas  for  the  A- 
coefficients  appearing  in  (79)  through  (82),  and  in  (83)  thro  igh  (86), 
are  the  following: 
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<■  j 

-  2  tt  p J r*bfh  F&  dr 
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1 
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i  : 

■  Zrrp J r*b  (j  -  /*.  dr  . 

The  terms  (79)  through  (82)  are  derived  from  the  generalized  forces  in 
ft.  and  <9  coordinates  due  to  the  quasi-steady  portion  and  the  apparent 

mass  contribution  (as  called  by  von  Karman  and  Sears  in  Reference  5) 

*  •  • 

of  the  aerodynamic  loads  resulting  from  hJ  ,  <9  and  6j  .  (The 
apparent  masses  associated  with  h-  and  04  are  negligible  compared 
with  the  corresponding  blade  inertial  properties.  )  In  other  words, 

(79)  and  (80)  are  the  rth  harmonic  cosine  and  sine  components, 
respectively,  of  the  generalized  force  in  />t  due  to  a  lift  distribution 


acting  at  quarter  chord  and  another  lift  distribution 

2-Tr/>b  n.r(b£_  Oj-Fq) 

'  j*i 

acting  at  half  chord,  while  (81)  and  (82)  are  those  in  &•  due  to  the 
same  two  lift  distributions.  (The  first  lift  distribution  is  entirely 
quasi-steady,  but  half  of  the  second  one  is  apparent  mass  contribution.  ) 


Based  on  airfoil  theory  of  two-dimensional  ‘hin  plate  and  assumed  T2r 
to  be  the  stream  velocity. 
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COMPUTER  PROGRAM 


Initially,  the*  computer  program  was  written  in  FORTRAN  IV  for  the 
IBM  704-1.  In  the  latter  stages  of  the  project,  the  program  was  simply 
converted  for  use  on  the  IBM  360/65;  no  effort  was  expended  to 
optimize  the  program  in  order  to  take  advantage  of  the  360/65  system 
capabilities.  Physically,  because  of  core  stoiage  limitations  of  the  7044 
and  for  convenience  in  solving  the-  equations  of  blade  motions,  the  program 
was  separated  into  six  parts.  Each  part  is  also  divided  into  several  sub¬ 
routines  for  further  convenience  in  programming,  checkout,  and  program 
modification.  The  major  steps  in  each  part  of  the  computer  program  are 
shown  in  the  flow  chart  presented  in  Figure  5.  Parts  1  and  2  of  the 
program  are  identical  to  those  presented  in  Reference  1,  since  the 
representations  of  the  blade  and  the  wake  wei  e  not  altered.  The  principal 
functions  of  each  part  of  the  program  are  listed  below. 

Part  1 

1.  Reads  the  necessary  inputs  for  Parts  1  and  2. 

2.  Computes  the  blade  and  wake  coordinates. 

Two  tape  drives  arc  used  in  this  part  of  the  program.  The  first 
is  used  for  temporary  storage  of  the  blade  and  wake  coordinates  used  in 
computing  the  induced  velocity  coeffic ient s,  and  the  second  is  used  to 
transfer  variables  and  the  matrix  of  induced  velocity  coefficients  which 
are  used  in  Part  2  of  the  program. 

Part  2 

1  Reads  the  necessary  inputs  for  Parts  2  and  3. 

2.  Computes  the  initial  trial  circulations  (represented  by  r  's 

n  Reference  1)  about  blade  segments  at  chosen  azimuth 
positions,  the  initial  blade  motions  and  other  variables 
used  in  Parts  2  and  3. 

3.  Normalizes  the  matrix  of  induced  velocity  coefficients. 

(The  normalized  matrix  is  used  in  solving  the  F  equations 

in  Part  3  .  ) 
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Three  tape  drives  are  used  in  this  part  of  the  program.  The  first 
tape  is  used  to  read  in  the  variables  and  the  matrix  of  induced  velocity 
coefficients  from  Part  1.  The  second  tape  stores  the  normalized  matrix 
of  induced  velocity  coefficients  used  in  solving  the  r  equations  in  Part  3; 
the  third  tape  transfers  the  initial  trial  r  ,  the  initial  trial  blade  motions 
and  other  variables  used  in  Part  3. 

Part  i 

1.  Solves  the  r  equations. 

2.  Solves  the  equations  of  blade  motions. 

3.  Computes  time  histories  of  blade  lift,  drag  and  pitching 
moment  distributions. 

4.  Computes  the  generalized  aerodynamic  forces. 

5.  Computes  the  time  histories  of  the  rotor  aerodynamic  forces 
and  moments  {due  to  one  blade)  used  in  Part  5. 

Four  tape  drives  are  used  in  Part  3  of  the  program.  The  first  is 
used  to  read  in  the  initial  trial  r  ,  the  initial  trial  blade  motions  and 
other  variables  used  in  Parts  3,  4,  and  5.  The  second  tape  retains  the 
normalized  matrix  of  induced  velocity  coefficients  used  in  solving  the 
equations.  The  third  tape,  which  has  been  generated  in  Part  6  of  the 
program,  is  used  to  retain  all  the  inverses  of  all  the  matrices  used  in 
solving  the  blade  equations  of  motions  for  each  harmonic  order.  The  fourth 
tape  is  used  to  store  the  time  histories  of  the  variables  to  be  harmonically 
analyzed  and  other  variables  used  in  Part  4  of  the  program.  The  time 
histories  of  the  rotor  aerodynamic  forces  and  moments  used  in  Part  5  of 
the  program  are  machine -punched  on  data  cards. 

Part  4 

Harmonically  analyzes  the  time  histories  of  the  variables 
generated  in  Part  3  of  the  program. 
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One  tape  drive  is  used  in  Part  4;  namely,  that  which  contains 
the  time  histories  of  the  variables  to  be  harmonically  analyzed. 

Part  5 

1.  Computes  the  harmonic  components  of  the  rotor  aerodynamic 
forces  and  moments  that  appear  in  the  trim  equations. 

2.  Computes  the  partial  derivatives  of  the  average  rotor  forces 
and  moments  with  respect  to  the  trim  constants.  (See(36).  ) 

3.  Computes  the  forces  and  moments  on  the  fuselage. 

4.  Computes  a  new  set  of  trim  variables. 

No  tape  drives  are  used  in  this  part  of  the  program.  The  time 
histories  of  the  rotor  aerodynamic  forces  and  moments  are  inputs  from 
data  cards  which  have  been  machine -punched  in  Part  3  of  the  program. 

The  integrals  used  in  Part  5  which  are  dependent  solely  on  the  blade 
geometric  and  inertial  properties  and  mode  shapes  are  computed  and 
machine-punched  in  Part  6  of  the  program  and,  thus,  inputed  on  data 
cards  in  Part  5. 

Part  6 

1.  Computes  the  coefficients  of  the  equations  of  blade  motions, 
(83)  through  (86),  to  form  a  matrix  for  each  harmonic  order. 

2.  Computes  the  inverses  of  the  matrices  which  are  used 
in  solving  the  equations  of  blade  motions  in  Part  3. 

3.  Computes  the  integrals  which  are  dependent  on  the  blades' 
geometric  and  inertial  properties  and  mode  shapes  to  be 
used  in  Parts  5  and  3. 

One  tape  drive  is  used  in  this  part  of  the  program;  namely,  that 
which  retains  the  inverses  of  the  matrices  and  the  integrals  appearing  in 
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Parts  1  and  6  oi  the  program  are  run  initially, and  the  binary  output 
tapes  from  each  part  are  saved.  Parts  1  and  3  are  then  run  sequentially 
a  total  of  five  times  to  provide  the  input  for  computing  the  rotor  fort  os 
and  moments  and  for  computing  the  partial  derivatives  of  the  rotor  forces 
and  moments  with  respect  to  each  of  the  trim  constants.  Part  S  is  then 
run  to  compute  a  new  set  of  trim  constants.  If  the  difference  between  the 
initial  set  of  trim  constants  and  the  new  set  is  not  within  prescribed 
limits,  Parts  2  and  3  are  then  rerun  with  the  new  set  of  trim  c  onstants 
to  provide  the  input  for  recomputing  the  rotor  forces  and  moments.  Part 
5  is  then  rerun  with  this  new  input  to  compute  a  third  set  of  trim  variables 
The  process  is  repeated  until  the  trim  variables  converge  within  prescribed 
limits.  When  the  process  has  converged,  the  output  tape  from  Part  3  is 
used  in  Part  4  to  harmonically  analyze  the  last-obtained  results. 

It  has  been  assumed  in  the  above  processes  that  the  partial 
derivatives  of  the  rotor  forces  and  moments  with  respect  to  the  trim 
constants  are  nearly  constant  within  the  range  under  consideration.  If 
this  were  not  true,  then  Part  3  would  need  to  be  rerun  to  obtain  a  new  set 
of  partial  derivatives  for  each  new  set  of  trim  constants. 
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38 


PART  4 


TAPE  FROM 
PART  3 


PART  5 


INTEGRALS 
FROM 
PART  6 


PART  6 


BLADE  GOEMETRIC, 
ELASTIC  A 
yMASS  PROPERTIES,, 
AND  VIBRATION 
MODES. 


HARMONIC 

ANALYSIS 


HARMONICS  OF 
LIFT.  DRAG 
A  PITCHING 
MOMENT 


I 


HARMONICS  OF 
rrt  AND 
GLAUERT 
COEFFICIENTS, 


FUSELAGE  MASS 
A  AERODYNAMIC 
PROPERTIES 


TIME  HISTORIES  OF 
iROTOR  AERODYNAMIC 
FORCES  A  MOMENTS  t 
FROM  PART  3 


INITIAL  TRIM 
CONSTANTS 


COMPUTE  PARTIAL 
DERIVATIVES  , 
DEFINED  BY  (36) 


COMPUTE  CERTAIN 
HARMONIC  COMPONENTS 
OF  ROTOR  AERODYNAMIC 
FORCES  A  MOMENTS 


COMPUTE  FUSELAGE 
FORCES  A 
MOMENTS 


COMPUTE  THE 
M-,  T-,  K-  AND 
2C-C0EFFICIENTS 


COMPUTE  THE 
A-COEFF I C I  ENTS 


COMPUTE  INTEGRALS 
APPEARING  IN  Cq. 

AND  INTEGRALS  ' 
USED  IN  PART  5 


TOW 

EACH  HARMONIC 
ORDER 


INVERT  THE 
MATRICES 


INPUTS 


COMPUTE  NEW 
SET  OF  TRIM 
CONSTANTS 


INPUTS 


NEW  SET  OF 
TRIM 

CONSTANTS 


I 


END 


INVERSES  OF 
.  MATRICES 
AND  INTEGRALS 
APPEARING 
IN  C. 

<1 


T 


OPTIONAL 

DATA 


I 


INTEGRALS 
USED  IN 
PART  5 


Figure  5  (CONTD).  COMPUTER  PROGRAM  FLOW  CHART. 
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COMPUTED  RESULTS  AND  COMPARISONS  WITH 


MEASURED  RESULTS 


The  computational  procedure  was  used  to  analyze  the  same 
four  flight  conditions  that  have  been  analyzed  in  Reference  1  They 
are  for  the  teetering  JH-1A  rotor  at  advance  ratios  0.26  and  0.08 
(identified  as  Flights  No.  67  and  I-Jo.  65,  respectively,  in  Reference  7) 
and  for  the  H-34  rotor,  which  consists  of  four  fully  articulated  blades, 
at  advance  ratios  0.29  (identified  as  Flight  No.  18  in  Reference  8)  and 
0.  18  (for  which  measured  data  are  obtained  similarly  to  those  i.. 
Reference  8  and  are  given  in  Table  III  of  Reference  13).  The  prolile 
drag  coefficients,  which  are  expressed  as  functions  of  angle  ol  attack 
and  Mach  number  in  the  computing  program,  are  based  on  data  from 
Reference  9  for  the  NACA  0015  airfoil  section  used  in  the  UH-1A  rotor 
blades  and  on  data  from  Reference  10  for  the  modified  NACA  0012  air¬ 
foil  section  used  in  the  H-34  rotor  blades.  The  fuselage  aerodynamic 
characteristics  are  obtained  from  Reference  11  for  UH-1A  and  Refer¬ 
ence  12  for  H - 34. 

The  azimuthal  variations  (or  the  time  histories  with  the  steals 
components  at  each  specified  radial  station  removed)  of  the  resulting 
blade  aerodynamic  loadings  (lift,  pitching  moment,  and  drag,  pe  r  unit 
span)  and  blade  bending  moments  are  plotted  for  a  number  ol  radial 
stations.  The  bending  moments  arc  obtained  by  summing  the 
contributions  from  the  responses  of  all  the  bending  modes  used.  They 
automatically  include  the  influence  of  inertia  loadings.  Available 
measured  data  (lift  distributions  and  flatwise  bending  moments)  are 
shown  in  the  appropriate  plots  lor  convenient  comparisons  in  Figures  6 
through  31.  For  the  lift  distribution  or  flatwise  bending  moment,  the 
radial  variation  of  the  steady  component  (or  time  average)  is  presented 
together  with  the  radial  variations  of  the  cosine  and  sine  components 
of  the  harmonics  in  blade  rotational  speed. 


Positive  flatwise  and  chordwise  bending  moments  tend  to  compress 
the  upper  surface  and  leading  edges,  respectively. 
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For  c.'u  h  flight  condition,  the  same  wake  model  and  the  same 
values  of  the  parameters  that  spec  lfy  the  wake  model  used  in  Refer¬ 
enda  I  are  used  here.  Five  flapwise  modes  (based  on  an  untwisted  blade 
with  zero  blade  angle)  are  used  in  Reference  1,  while  live  or  six  bending 
modes  (in  which  both  flapwise  and  chordwise  deflections  are  admitted 
and  twist  effects  are  included)  and  two  torsional  modes  art*  used  here. 

The  computed  lift  and  flapwise  bending  moment  results  obtained  here 
are  not  significantly  closer  to  the  measured  results  than  those  obtained 
in  Reference  1.  Therefore,  the  establishment  of  both  a  better  wake 
model  and  further  refinements  in  the  prediction  of  blade  motions  may 
be  necessary  for  the  configurations  studied.  However,  the  present 
results  do  demonstrate  that  the  computational  procedure  developed 
by  the  present  eftort  is  suitable  for  taking  additional  blade  motions  into 
acc  ount  and  can  be  used  to  yield  trim  constants. 

Reasonable  agreement  between  measured  and  computed  lift 
distributions  is  obtained  (Figures  6,  7,  13,  14,  19,  20,  26,  27), 
but  the  flatwise  bending  moment  agreement  is  less  satisfactory  (Figures 
8,  d,  15,  16,  21,  22,  28,  29).  The  disagreement  in  the  flatwise 
bending  moment  values  may  be  the  result  of  an  accuracy  problem.  The 
moments  are  relatively  small, and  they  represent  small  differences  of 
large  numbers  (lift  moments  and  inertial  force  moments)  for  the  con¬ 
figurations  studied.  In  all  cases,  both  the  measured  and  computed  llat- 
wise  bending  moments  on  any  blade  section  have  magnitudes  less  than  the 
moment  acting  on  some  inboard  blade  section  due  to  the  blade  weight 
alone,  while  the  lift  distribution  over  the  outboard  portion  of  the  blade 
goes  as  high  as  20  to  30  times  the  blade  weight  distribution.  Thus,  small 
differences  between  the  measured  and  computed  lift  distributions  can 
result  in  large  percentage  differences  in  flatwise  bending  moments. 
Similarly,  small  errors  between  computed  and  actual  natural  frequencies 
lead  to  large  errors  in  inertial  forces  at  forcing  frequencies  near  the 
resonance  frequencies. 

Five  bending  modes  are  used  in  the  UH-1A  cases  and  six  in  the  H-34 
casts,  with  the  natural  frequency  of  the  highest  order  mode  slightly 
over  ten  times  the  rotor  shaft  rotational  speed.  (The  present  computer 
program  can  handle  a  total  of  ten  bending  and  torsional  modes.  ) 
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Calculated  drag  distributions  and  pitching  moment  distributions 
are  shown  on  Figures  10,  11,  17,  18,  23,  24 ,  30,  and  3 1 . 

Chordwise  bending  moments  are  shown  on  Figures  12  (UH-1A) 
and  25  (H-34).  The  chordwise  bending  moments  as  computed  by  summing 
the  contributions  from  the  responses  of  all  the  bending  modes  used  do 
not  agree  with  estimates  based  on  computed  blade  drag  loading.  Super¬ 
position  of  natural  modes  is  ineffective  for  the  cases  considered  because 
the  blade  chordwise  elastic  deflection  is  small  and  is  almost  all  due  to 
a  /ingle  bending  mode  which  is  dominated  by  chordwise  deflection. 

Thus,  in  effect,  the  chordwise  deflection  distribution  is  approximated 
by  only  one  curve  which  can  have  a  curvature  (or  bending  moment)  dis¬ 
tribution  quite  different  from  the  actual  one.  That  is,  unless  the  loading 
has  the  same  shape  as  the  inertial  loading  in  that  single  natural  mode, 
the  corresponding  moment  distributions  will  not  have  the  same  shape. 

Measured  chordwise  bending  moments  are  available  at  four 
stations  for  the  H-34  case.  Both  the  measured  and  computed  chord- 
wise  bending  moments  for  the  H-34  rotor  at  an  advance  ratio  of  0.  29 
are  shown  in  Figure  25.  The  poor  agreement  is  attributed  to  the  factors 
given  above,  although  it  must  be  admitted  that  more  cases  should  be 
investigated  before  a  strong  conclusion  can  be  formed.  By  considering 
that  the  moment  about  the  lead-lag  hinge  due  to  chordwise  centrifugal 
loading  cancels  that  due  to  the  computed  drag,  the  resulting  steady 
inplane  bending  moment  distributions  have  been  computed  by  hand.  They  are 
shown  in  the  table  on  page  43,  together  with  the  measured  and  computed 
steady  chordwise  bending  moments,  all  in  inch-pounds.  The  measured 
value  of  - 1  392 1  at  r/R  -  0.  825  corresponds  to  a  uniform  drag  loading  on 
the  blade  outboard  of  that  station  of  about  97  pounds  per  foot.,  which  can 
hardly  be  considered  reasonable.  The  disagreement  between  the  two  sets 
of  computed  results  is  disturbing  (this  would  be  true  even  if  no  measured 
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All  the  other  natural  bending  modes  are  dominated  by  the  flatwise  motion. 
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STEADY 

TABLE 

INPLANE  BENDING  MOMENTS,  H- 

34,  0.29 

r/R 

STEADY  INPLANE 
BENDING  MOMENT 
COMPUTED  FROM 
BLADE  LOADING 

MEASURED 

STEADY  CH0RDWISE 
BENDING  MOMENT 

COMPUTED 

STEADY  CH0RDWISE 
BENDING  MOMENT 

0.150 

-2670 

2379 

390 

0.275 

-4530 

mm 

291 

0.375 

-4600 

3776 

190 

0.575 

-3293 

-4416 

-5 

0.125 

-620 

-13921 

IS 

values  wei  available).  In  effect,  the  calculated  results  suggest 
that  it  would  be  advisable  to  calculate  steady  deflections  (and  moments) 
from  the  steady  loads  and  to  use  the  modal  description  to  account  for 
the  time  varying  deflections  (and  moments). 

The  calculated  and  measured  trim  data  are  given  below  for  each 
flight  case  considered.  It  should  be  recalled  that  the  gross  weight,  c.  g. 
position,  fuselage  aerodynamic  forces  (lift,  drag,  pitching  moment), 
rotor  rotational  speed,  and  flight  velocity  are  input  data.  Direct  com¬ 
parison  of  the  angle  components  (a,  )  calculated  with  those 

measured  should  be  made  with  caution  for  several  reasons.  First, 
aerodynamic  effects  usually  depend  on  a  combination  of  the  components; 
e.  g.  ,  location  of  the  "control  plane".  Second,  the  accuracy  of  the  experi¬ 
mental  data  has  not  been  specified.  Third,  the  induced  angles  of  attack 
were  not  measured  and  they  can  have  magnitudes  comparable  to  the 
magnitudes  of  some  of  the  control  angles.  Fourth,  the  structural 
deformation  effects  have  not  been  explicitly  shown,  and  these  can  be 
important. 
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UH-1A  at  ^  =  0,  26 


The  computed  trim  constants  in  degrees  art*  listed  below  (with  the 
corresponding  measured  values,  obtained  from  Reference  7,  shown  in 
parentheses ): 


<*s  = 

5.  35 

deg 

(  6.  50, 

measured) 

^  = 

19.  03 

deg 

(19.  20, 

measured) 

&1C  - 

1.  10 

deg 

(  2. 86, 

measured) 

Sts  = 

-6.  47 

deg 

(-7.  34, 

measured) 

The  computed  average  rotor  forces  are 

thrust  =  ~  Z  =  6,610  pounds 

longitudinal  force  =  X  =  -535  pounds 

— 

side  force  =  ~  Y  =  -271  pounds 

yj  _ 

torque  =  @  =7,250  foot-pounds 

rtf, 

pitching  moment  =  -  606  foot-pounds 

In  the  expression  of  steady  plus  cyclic  blade  flapping  responses, 

A  *  fio  +  fitc  cos  Sin 

The  results,  all  given  in  degrees,  are 

6  -  0  deg  (-0.07,  measured) 

3 /c  -  0.24  deg  (-1.10,  measured) 

3,  =-0.23  deg  (  2.39,  measured) 
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UH-1A  at  M  -  °-  08 

U  s  i  ng  the  same  notations  and  units  as  in  the  previous  case,  the 
results  for  this  case  are 


OCs  = 
$o  ~ 
$tc  - 
&ts  -- 


5.  03  deg 
1 7.  57  deg 
0.81  deg 
-2.  56  deg 


(  4.  SO,  measured) 
(14.97,  measured) 
(  1.49,  measured) 
(-  1 .  73 ,  measured) 


thrust  =  6,620  pounds 
longitudinal  force  =  -396  pounds 

side  force  =  -203  pounds 

torque  =  5,  420  foot-pounds 
pitching  moment  =  3S8  foot-pounds 


0  deg 
0.83  deg 
-l.  53  deg 


(  0.  16,  measured) 
(-2.07,  measured) 
( -0.  33,  measured) 


H-  34  at  /x  =  0.  29 

The  following  results  are  obtained  for  the  present  case: 


a*  = 

5.  29  deg 

(  5.40, 

measured) 

= 

1 6. 04  deg 

(17.  35, 

measured) 

*tc  = 

2.24  deg 

(  2.61, 

measured) 

= 

-6. 37  deg 

(-7. 82, 

mea  sured) 

thrust  = 
longitudinal  force  = 
side  force  = 
torque  = 
pitching  moment  = 


11,630  pounds 
46  pounds 
-721  pounds 
23,780  foot-pounds 
372  foot-pounds 
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Ac 

- 

3.  82 

deg 

(  3.71, 

measured) 

- 

-0.  01 

deg 

(-0. 02, 

measured) 

firs 

- 

o 

00 

deg 

(-0. 24, 

measured) 

- 

8.  16 

deg 

(  9.  15, 

measured) 

he 

- 

0.  16 

dog 

(-0. 31, 

measured) 

ZfS 

- 

0.  10 

deg 

(  0.  14, 

measured) 

The  last  three  items,  all  in  degrees,  are  the  quantities  as  defined  in  the 
expression  of  steady  ulus  cyclic  blade  lagging  responses; 

A  ■  *  Tfc  cos  *■  T rs 


H-34  atM  =  0.  18 

(  2.  50,  measured) 
(13.  69,  measured) 
(  1.19,  measured) 
(-4.  70,  measured) 


ors  =  2.  37  deg 

e0  =  13.  12  deg 

Gk  =  3.  05  deg 

9,s  -  -3.  25  deg 


thrust  = 
longitudinal  force  = 
side  force  = 
torque  = 
pitching  moment  = 


1  1,820  pounds 
81  pounds 
-444  pounds 
14,640  foot-pounds 
646  foot-pounds 


A  = 

3. 

86 

deg 

(  3. 

93, 

measured) 

Arc  ~ 

-0. 

07 

deg 

(  o. 

61, 

mea  sur  ed) 

A/s  = 

0. 

55 

deg 

(-0. 

52, 

measured) 

-K  — 

Lo 

5. 

01 

deg 

(  5- 

99, 

mea  sured) 

he  = 

0. 

00 

deg 

(-0. 

13, 

measured) 

hs  - 

0. 

1 1 

deg 

(  o. 

04, 

mea  sured) 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A  computational  procedure  has  been  developed  that  predicts  rotor- 
blade  aerodynamic  loadings,  bending  moments,  and  realistic  blade  motions, 
as  well  as  trim  constants.  The  scope  of  the  blade  motion  analysis  of 
Reference  1  was  extended  by  incorporation  of  lead-lag  and  edgewise  bend¬ 
ing  and  torsional  degrees  of  freedom.  Although  these  additional  blade 
motion  degrees  of  freedom  have  been  included,  the  computed  lift  and  flap- 
wise  bending  moment  distributions  for  the  UH-1A  and  CH-34  rotors 
obtained  here  are  not  significantly  closer  to  the  measured  results  than 
those  obtained  in  Reference  1.  Therefore,  it  appears  that,  at  least  for 
these  configurations,  further  improvement  in  the  correlations  will  probably 
require  the  establishment  of  a  better  wake  model  and  an  even  more  complete 
representation  of  the  blade  motions. 

In  both  Reference  1  and  the  present  work,  a  fixed  skewed  helix  was 
assumed  for  the  location  of  the  vorticity  in  the  wake.  It  appears  that  a 
realistic  wake  model  is  needed  and,  to  this  end,  both  theoretical  and 
experimental  studies  of  the  disposition  of  vortical  elements  in  rotor  wakes 
are  recommended.  An  experimental  program  alone  probably  will  not 
suffice  since  quantitative  data  are  required  which  can  be  extended  to  any 
desired  rotor  configuration.  Current  experimental  techniques  are 
susceptible  to  errors  which  are  not  adequately  predictable.  For  instance, 
the  presence  of  boundaries  (such  as  wind  tunnel  walls)  must  be  accounted 
for  in  the  interpretation  of  experimental  data.  Direct  application  of  the 
wake  positions  measured  on  one  configuration  might  introduce  serious 
errors.  A  concommitant  theoretical  effort  is  recommended,  therefore, 
to  enable  evaluation  of  experimental  techniques  and  to  develop  a  means 
of  theoretically  predicting  the  wake  location  for  any  specified  rotor  system. 

Further  improvement  of  the  correlations,  particularly  in  the 
bending  moments,  might  necessitate  a  better  representation  of  the  blade 
motions.  Configurations  such  as  the  UH-1A  and  CH-34  rotors  are  designed 
to  maintain  low  bending  moments  of  the  blades.  Consequently,  smill 
errors  in  predicting  the  blade  motions  and  the  resultant  lift  can  produce 
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relatively  large  errors  in  the  computed  bending  moments,  good  prediction 
of  the  bending  moments  might  depend,  therefore,  on  a  precise  description 
of  the  inertial  and  elastic  characteristics  of  the  blade  as  well  as  the 
definition  of  the  wake. 

Finally,  it  is  recommended  that,  if  this  computer  program  is  to 
see  practical  utilization  it  should  be  revised  to  make  efficient  use  of  the 
capabilities  of  the  IBM  360/65  on  which  it  is  now  being  run.  This  computer 
program  was  written  initially  in  FORTRAN  IV  for  the  7044.  Later  in  the 
project,  the  program  was  simply  converted  to  the  IBM  360/65.  There  is 
no  effort  to  optimize  the  program  to  the  new  equipment. 

In  summary,  the  following  further  efforts  are  recommended:  (a) 
investigate  means  for  improving  the  wake  model,  (b)  determine  the 
relative  advantages  of  alternative  representations  of  blade  motion 
equations,  and  (c)  refine  the  computer  program  in  order  to  take 
advantage  of  the  IBM  360/65  system  capabilities. 
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LIFT  -  lb/ft 


400 t  r/R  =  0.95 

STEADY  COMPONENT 

=  288  Ib/ft,  MEASURED 
=  294  lb/ft,  COMPUTED 


£  o 


AZIMUTH  ANGLE  -  degrees 


Figure  6.  MEASURED  AND  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
LIFT  DISTRIBUTION;  UH-1A  AT /4  =  0.26. 
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LIFT  -  lb/ft  LIFT  -  Ib/ft 


STEADY  COMPONENT 


-  COMPUTED 

0  MEASURED 


COSINE  COMPONENT  SINE  COMPONENT 
1st  HARMONIC 


MEASURED  AND  COMPUTED  HARMONICS  OF  BLADE 
LIFT  DISTRIBUTION;  UH-1A  AT  p  =  0.26. 


BENDING  MOMENT  -  in. -It 


0  90  180  270  360  0  90  180  270  360 


AZIMUTH  ANGLE  -  degrees  AZIMUTH  ANGLE  -  degrees 

Figure  8.  MEASURED  AND  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
FLATWISE  BENDING  MOMENTS;  UH-1 A  AT  px  =  0.26. 
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Figure  9.  MEASURED  AND  COMPUTED  HARMONICS  OF  BLADE  FLATWISE 
BENDING  MOMENT;  UH-IA  AT/c  =  0.26. 


r/R  =  0.40 

STEADY  COMPONENT  =  10.5  Ib/ft 


i - - - r . . . i  . i 


r/R  =  0.85 

STEADY  COMPONENT  =  16.8  lb/ft 


Figure  10.  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
DRAG  DISTRIBUTION;  UH-1A  AT /x.  =  0.26. 


PITCHING  MOMENT  -  ft-lb/ft 


Figure  II.  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE  PITCHING  MOMENT 
(ABOUT  MIDCHORD)  DISTRIBUTION;  UH-IA  AT =  0.26. 
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STEADY  COMPONENT 
=  6640  in. -lb 


-2 ' 


i  r 
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r/R  =  0.80 


STEADY  COMPONENT 
=  1150  in. -lb 


1  r 


r/R  =  0.925 


STEADY  COMPONENT 

=  210  in. -lb 


90  180  270 

AZIMUTH  ANGLE  -  degree* 


360 


Figure  12.  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE  CHORDWI SE 
BENDING  MOMENTS;  UH-1A  AT/x  =  0.26. 
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LIFT  -Ib/ft 


r/R  =  0.90 


AZIMUTH  ANGLE  -  degrees 


Figure  13.  MEASURED  AND  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
LIFT  DISTRIBUTION;  UH-JA  AT  /*  =  0.08. 
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STEADY  COMPONENT  COSINE  COMPONENT  SINE  COMPONENT 


r/R  r/R  r/R  r/R 

Figure  14.  MEASURED  AND  COMPUTED  HARMON KS  OF  BLADE 
LIFT  DISTRIBUTION;  UH-1A  AT  /* =  0.08. 
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=  -10149  i n. -1 b, 'MEASURED 


y 

STEADY  COMPONENT 

=  -420  in. -1b,  MEASURED 
=  -221  in. -lb,  COMPUTED 


=  -140  in. -1b.  MEASURED 


=  -307  in. -1b,  MEASURED 
=  -401  in. -1b,  COMPUTED 


AZIMUTH  ANGLE  -  degrees 


-  MEASURED 

L - COMPUTED 

S  EADY  COMPONENT 

=  -2831  in. -1b.  MEASURED 
=  _635  in. -lb.  COMPUTED 


=  -3218  in. -1b.  MEASURED 
=  -680  in. -lb,  COMPUTED 


=  -1601  in. -lb.  MEASURED 
=  -174  in. -1b,  COMPUTED 


AZIMUTH  ANGLE  -  degrees 


Figure  15.  MEASURED  AND  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
FLATWISE  BENDING  MOMENTS:  UH-IA  AT  M  =  0.08. 
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BENDING  MOMENT;  UH-1A  AT/*  =  0.08. 
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DRAG  -  lb/ft 


0  90  180  270  360 


AZIMUTH  ANGLE  -  degree* 


0  90  ISO  270  360 


AZIMUTH  ANGLE  -  degree* 


Figure  17.  COHPUTEO  AZIMUTHAL  VARIATIONS  OF  BLADE 
DRAG  DISTRIBUTION;  UH-IA  AT =  0.08. 
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Figure  18.  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE  PITCHING  MOMENT 
(ABOUT  MIDCHORD)  DISTRIBUTION;  UH-IA  AT  ^  =  0.08. 
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=  180  Ib/ft.  MEASURED 
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=  89  Ib/ft,  MEASURED 
=  54  Ib/ft,  COMPUTED 
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f=  245  Ib/ft.  MEASURED 
=  225  Ib/ft,  COMPUTED 
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\  '▼s  232  lb/ft,  MEASURED 
W  =  237  Ib/ft,  COMPUTED 


90  180  270  360  0  90  180  270  360 


AZIMUTH  ANGLE  -  degree* 


AZIMUTH  ANGLE  -  degrees 


r/R  =  0.95 


C  0 


£  -100-  I  1/ 

3  V.  //  STEADY  COMPONENT 

A/  =  229  Ib/ft,  MEASURED 

=  241  Ib/ft,  COMPUTED 

-200*  r  1  ■  ■  >  t 

90  180  270  360 

AZIMUTH  ANGLE  -  degrees 

Figure  19.  MEASURED  AND  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
LIFT  DISTRIBUTION;  H-34  AT  p  =  0.29. 
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LIFT  -  Ib/ft 


Figure  20.  MEASURED  AND  COMPUTED  HARMONICS  OF  BLADE 
LIFT  DISTRIBUTION;  H-34  AT,u  =  0.29. 
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BENDING  MOMENT 


0  90  180  270  360  0  90  180  270  360 

AZIMUTH  ANGLE  -  degrees  AZIMUTH  ANGLE  -  degrees 


Figure  21.  MEASURED  AND  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
FLATWISE  BENDING  MOMENTS;  H-34  AT//  =  0.29. 
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r* -  t  1  ■  t  i - » 

0  90  180  270  360 

AZIMUTH  ANGLE  -  degrees 


Figure  23.  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
DRAG  DISTRIBUTION;  H-34  AT  ju.  =  0.29. 
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0  90  180  270  360 

AZIMUTH  ANGLE  -  degrees 


Figure  24.  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE  PITCHING  MOMENT 
(ABOUT  MIDCHORD)  DISTRIBUTION;  H-34  AT // =  0.29. 
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BENDING  MOMENT 


AZIMUTH  ANGLE  -  degrees  AZIMUTH  ANGLE  -  degrees 


AZIMUTH  ANGLE  -  degrees 


Figure  25.  MEASURED  AND  COMPUTED  AZIMUTHAL  VARIATIONS  OF 

BLADE  CHORDWISE  BENDING  MOMENTS;  H-3H  AT/4  =  0.29. 
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Figure  26.  MEASURED  AND  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
LIFT  DISTRIBUTION:  H-34  at  H  =  0.18. 
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Figure  27.  MEASURED  AND  COMPUTED  HARMONICS  OF  BLADE 
Ufr  DISTRIBUTION;  H-34  AT  fu  =  0.18. 
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Figure  28.  MEASURED  AND  COMPUTED  AZIMUTHAL  VARIATION  OF  BLADE 
FLATWISE  BENDING  MOMENTS;  H-34  AT  =  0.18. 
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BENDING  MOMENT  -  in. -lb  *  10”a  BENDING  MOMENT 


STEADY  COMPONENT  COSINE  COMPONENT  SINE  COMPONENT 


r/«  r/»  f/N  r/R 


Figure  29.  MEASURED  AND  COMPUTED  HARMONICS  OF  BLADE  FLATWISE 
BENDING  MOMENT;  H-34  AT /.<  =  0.18. 
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Figure  30.  COMPUTED  AZIMUTHAL  VARIATIONS  OF  BLADE 
DRAG  DISTRIBUTION;  H-34  AT jul  =  0.18. 
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AZIMUTH  ANGLE  -  decrees 

Figure  31.  COMPUTEO  AZIMUTHAL  VARIATIONS  OF  BLADE  PITCHING  MOMENT 
(ABOUT  MIDCHORD)  DISTRIBUTION;  H-34  AT/*  =  0.18. 
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APPENDIX 


ORTHOGONALITY  RELATIONS  OF  BENDING  VIBRATION  MODES 
OF  A  TWISTED  ROTATING  BEAM 


The  bending  deformation  of  a  twisted  (rotating  or  nonrotating) 
beam  is  generally  not  in  one  plane.  For  this  reason,  recent  practice 
in  load  and  deformation  analysts  of  a  twisted  propeller  or  helicopter 
blade  utilizes  bending  vibration  modes  of  the  twisted  rotating  blade 
with  both  flatwise  and  c’  ordwise  deflection  components  as  generalized 
coordinates.  Since  these  modes  all  satisfy  the  same  two  governing 
differential  equations  and  the  same  root  and  tip  conditions,  they  satisfy 
certain  relations  (usually  called  the  orthogonality  relations)  among  them 
selves.  This  note  presents  a  self-contained  derivation  of  the  general 
equations  f-om  which  the  relations  pertaining  to  blades  with  given  root 
conditions  can  be  obtained  immediately  upon  using  the  root  conditions. 

1  he  relations  are  useful  in  detecting  the  accuracy  of  results  obtained 
irom  vibration  analysis  and  in  simplifying  the  equations  of  blade 
motions  in  which  certain  bending  vibration  modes  are  used  as  a  group 
of  the  generalized  coordinates.  A  blade  with  flapping  and  lead-lag  hinges 
is  among  those  used  in  the  examples. 

Existing  methods  for  bending  vibration  analysis  of  twisted 
rotating  blades  (such  as  Reference  4)  assume  that  the  undeformecl 
blade  elastic  axis  is  a  straight  line  which  lies  in  the  plane  of  rotation 
and  intersects  the  axis  of  rotation.  Furthermore,  the  c.  g.  and  central 
polar  axes  arc  assumed  to  be  coincident  with  the  elastic  axis,  and  the 
built-in  twist  is  assumed  to  be  about  the  elastic  axis.  All  these 
assumptions  will  be  made  here.  In  blade  load  and  deformation  analysis, 
Lagrange's  equations  of  blade  motions  take  the  actual  axis  offsets 
and  built-in  loning  into  a  .count  by  having  constant  generalized  force 
terms  that  depend  on  the  offsets  and  built-in  coning.  fThe  offset 

7  7 


between  the  elastic  and  c.  g.  axes  also  causes  coupling  between  bending 
and  torsional  modes.  ) 

Consider  a  blade  section  at  a  distance  r  from  the  axis  of  rotation, 


> , 

C3n 


about  which  the  blade  (as  well  as  the  system)  is  rotating  in 

vacuum  at  a  constant  angular  velocity  Jl  .  With  reference  to  Figure  3<i. 
the  j.-axis  is  the  axis  ol  rotation,  the  r-axis  is  perpendicular  to  the 
paper  and  pointing  to  the  reader  and  the  r  -axis  is  perpendicular  to 
the  ^ y  -plane.  When  the  blade  is  undeformed,  the  section  is  shown  by 
the  solid  line:  its  elastic  center  lies  on  the  r  -axis  and  its  major 
principal  central  xis  makes  an  angle  0(r  ),  called  the  local  blade 
ang  •,  with  the  plane  of  rotation  (i.  e.  ,  the  -plane).  With  the 
presence  ot  b^r  ing  deformation,  the  same  section  is  shown  by  the 


78 


/ 


dotted  line:  its  elastic  center  is  displaced  by  (  «  ,  is-  )  or  (  c  ,  h  )  in  a 
plane  parallel  to  the  q  y  -plane,  while  6  is  assumed  to  remain  unchanged. 
The  deflection  components  u  (  r  )  and  v'(r)  are  called  the  chordwise 
and  flatwise  bending  deflection  components,  respectively,  while  c(r) 
and/i(r)  are  called  the  inplane  and  flapwise  components,  respectively. 
From  Figure  3Z, 


or 


h  =  VCOS  Q  +  lC.  Sir)  0 
c  ■  u  cos  6  -  ir  str)  6 


ir  = 

h 

cos 

O 

-  C  Sir 7 

u.  - 

c 

cos 

0 

■h  h  Sin 

(87) 


(88) 


Assuming  that  the  blade  sections  remain  plane  during  bending, 
the  small  deflection  described  above  is  the  result  of  the  small  rotation 
of  each  blade  section  about  an  axis  which  lies  in  the  plane  of  the 
section  and  intersects  the  elastic  axis.  The  rotation  has  one  component 


dh 


equal  to  radians  in  the  negative  ^-direction  and  the  other  component 
equal  to  dJdL  radians  in  the  ^.-direction,  according  to  the  right-hand 
rule  of  representing  rotations  by  vectors.  The  curvature  of  the 


projection  of  the  deformed  elastic  axis  on  the  ^r-plane  or  on  the  rrj- 

plane  is  equal  to  -y-j  or  %  respectively. 

dr  dr 


Let  G-  or  F  bf  the  curvature  of  the  projection  of  the  elastic  axis 
on  the  plane  that  contains  the  local  major  principal  central  axis  and 
the  r  -axis  or  on  the  plane  that  contains  the  local  minor  principal 
central  axis  and  the  r-axis,  respectively  Also,  let  M  and  W  be  the 
flatwise  (about  the  local  major  principal  central  axis)  and  chordwise 
(about  the  local  minor  principal  central  axis)  bending  moments  acting 
on  a  blade  section.  According  to  the  engineering  beam  theory,  which 
assumes  that  the  beam  cross  sections  remain  plane  during  bending,  the 


Positive  M  and  M  tend  to  compress  the  upper  surface  and  leading 
edge,  respectively. 
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following  relations  hold  at  any  blade  section: 


FI  =  8  F 


(89) 


and 


vl 


b  G 


(90) 


where  8  and  t>  ,  with  b>5  usually,  are  the  flatwise  and  chordw'ise 
flexural  rigidities  of  the  blade  section. 

On  account  of  the  blade's  being  twisted  (namely,  ~  f  0),  the 
curvatures  F  and  6  (each  belongs  to  the  projection  of  the  elastic  axis 
on  a  plane  that  rotates,  as  r  changes,  about  the  r  -axis  with  respect 
to  a  coordinate  system  fixed  to  the  blade)  are  not  given  by 


d*  V- 
drl 


and 


,2 

a  u 
dr * 


but  are  given  by 


d*h  -  d *  c  „ 

— — j  cos  &  ~  - “r  -sto  6 

dr  dr 

J  c  cos  e  *-  ^4 


dr 


d  r ' 


or,  by  the  use  of  (87), 


<?  = 


d*ir 

de 

r  6  - 

dec 

d'O 

f-  - 

dr* 

Jtr 

dr 

dr* 

\  dr) 

d*u 

,  dO 

dt~ 

dZF 

sdPd- 
v-  -  -  \  a 

d  r 

dr 

d  r 

d  r  * 

K  dr 

(01) 


(9<2) 


Let  A/ and  M  be  the  flapwise  and  inplane  components  of  the 
bending  moment  acting  on  a  blade  section  such  that 
V7  -  M  cos  0  r  M  sen  G 

-  &F  cos  6  i*-  bCrStnO 

-  3  ~T~l  cos  0  — —  sen  &\  cos  6 

'  dr  dr 

■b  b{~£~2  COS  0  sen  Q 


d  r ' 


(03) 
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and 


M  M  cos  &  v?  sen  0 


-  0  G  COi  C?  Br  Sen 


etc  -  d'h 

0  (  i  COS  b>  — 2  5,;/?  f  1  b 

V  dr  dr 


B  U  h  d  C 

-  5  f— ;  *  ^6?5  <9 - -  $m  &)A£n  H  . 

dr  d  > 


(94) 


Let  4;  and  p  be  the  components  of  loading,  perpendicular  to  the 
deformed  elastic  axis,  per  unit  length  of  blade  in  the  j  -  and  p- 
directions,  respectively,  Then,  for  a  rotating  blade  in  free  bending 
oscillations  at  a  natural  frequency  of  00  radians /  second  ,  the  loading 
component  p  is  given  by 


where  u  is  the  mass  per  unit  length  of  blade  and,  with  R.  being  the 
blade  tip  radius, 


/  - 


(96) 


is  the  centrifugal  force  acting  on  a  blade  section  at  any  r  . 
the  centrifugal  force  has  a  shear  component  equal  to 


X 


d  h 
aLr 


Since 


the  upward  component  of  transverse  loading,  or  the  contribution  to 
o  ,  due  to  CL  is  equal  to  the  last  term  of  (°5).  Carrying  out  the 


The  Coriolis  force,  which  has  the  effect  of  increasing  A  by 

-  dfl  J  C  u  dr  , 

is  neglected  in  vibration  analysis. 
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differentiation,  (95)  becomes 


J*h 


_  y ,  cX  n  2  Uh 

p  ~  on  h  t~  r  ~ — ~t  u  n  r  — —  . 

dr  a.  r 


(97) 


Similarly,  the  corresponding  p  can  be  expressed  as 


-  2  d  C  t  dc 

p  =  u  u>  c  *  d  -  p  Q.  r  — —  *  p  Cl  c 
ct  r 


(98) 


The  last  term  of  (98)  is  due  to  the  fact  that  when  the  c.  g.  of  a  mass 
is  off  from  the  r  -axis  in  the  y  -direction,  there  is  a  centrifugal  force 
component  in  the  y  -direction  acting  on  the  mass.  The  governing 
equations  for  free  vibrations  are 


d  V 

TP 


‘  p 


(9Q) 


and 

T ' M 

tt  * p  (,00) 

with  M ,  Pd  ,  p  and  p  as  given  by  (93),  (94),  (97),  and  (98),  respectively. 
When  any  particular  cth  or  jth  mode  is  under  consideration,  the 
subscript  t  or  j  will  be  added  to  uJ,  h  ,  c  ,  p  ,  p  ,  M  ,  f  and  0  • 

The  four  boundary  conditions  at  the  blade  tip,  where  r  =  R  , 


(99)  and  (100)  are  in  agreement  with  Reference  6,  since  the  two 
scalar  equations  below  (4.  8)  of  Reference  6  are  the  following  two 
combinations  of  (99)  and  (100): 

(99)coS  0  -(1 00)^0 

and 

(100)  cos0  a  (99)  . 
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may  be  written  as 


(")r  * 


(101) 


Since  (99)  and  (100)  are  a  pair  of  fourth-order  differential  equations  in 
h  and  c  ,  four  more  boundary  conditions  are  required  for  each 
actual  case.  These  additional  boundary  conditions  are  the  "root" 
conditions  which,  for  the  present  analysis,  are  assumed  to  be 
specified  at  r  =  r0  . 


By  first  multiplying  (99)  and  (100)  for  any  mode  by  h  and  C 
of  any  mode,  respectively,  then  adding  them  together,  and  finally 
integrating  over  *'  from  f0  to  R  , 

{ 

r* 

V  (pt  hj  +  pc.)  dr  .  (1 0?. ) 

O 

By  integrating  by  parts  twice,  using  ( 1  0  1 ),  the  left-hand  side  of  (102) 
may  be  expressed  as 

/  (  *  h,  + — -r-T-  C,  }dr 

Or  (  dr1  J  dr*  J  ) 


By  the  use  of  (96)  through  (98f, 

J  (P^j  +PiCj)dr 

ro 

*  “>*f  MWj  +c.cj)dr 

+ 1  {£(x&fa*£(x7?)c4dr 


(103) 


+  a 


7  ^c^r, 
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which,  by  applying  integration  by  parts  to  the  integral  that  contains  \ 
and  subsequently  using  (96),  may  be  written 


ZuJiJ  )  dr  +  n1  jj  juC^jdr 

\  e 

'  ijl  *rdr/ '(%  It  '  27  it)*’- +1dCX.r,  {  MrtJr\ 

Now,  substitution  of  (103)  and  (104)  into  (102)  yields 


(104) 


Kcj  ~nXtj  -\jJr  Mrdr)  +  3(J  =  u>* 

r  O 


(105) 


where 


M,  d  cl- 1  dr 


*/M 


M^)dr 


) 


*  /  *  ~  ~ 


K  fy)dr 


(106) 


A  A* 

rcj  £  J 

ro 

-/,;a>")(£  &  -  %- 

-d^dr 

dr 

(107) 

K  *J*a  M,  * c< cj) dr  '  JX u(K  ^ 

ro  0 

ru.L  Uj)dr 

(108) 

.  dh,  dc  ) 

A‘j  drCir.r, 

(109) 

B  =/m  »,  +at  dc< 

lJ  1  *  air  dr  J  1  dr  dr  J  J 

r-rc 

(HO) 

Since  ( 1  02 ) and.hencc,  ( 1  05)  hold  when  c  and  j  are  interchanged  and 
since 


jt 


lJ 


r^ 


A 

T 

V 


and  Mjg  -  MtJ  , 
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it  is  clear  that 


/ 


AT  •  1  : 


s  r  a 

'*  /L  r 


ir  t-  3,  M 


J  ‘■J 


(111) 


Subtracting  (111)  from  ( 1  05), 


(,c2  -  UL.?  M- ■  [5  & 

C  J  '  tj  OJ  „ 


U  -  *<*„  W' 


•_r  dr  -  O 


(112) 


Equations  (1 1  2 )  and  (105)  are  the  desired  equations.  Some  root 
constraints  of  practical  interest  will  be  dealt  with  in  the  examples 
that  follow.  As  the  examples  show,  (111)  yields  a  ’•elation  stating  that 
there  is  no  inertial  coupling  between  any  two  different  modes. 
Consequently,  (105)  yields  one  relation  stating  that  the  elastic  and 
centrifugal  couplings  between  any  two  different  modes  cancel  each 
other,  and  another  relation  stating  that  for  each  mode,  the  strain 
energy  at  either  extreme  configuration  minus  the  work  done  by  the 
centrifugal  loading  during  the  change  of  configuration  from  the  neutral 
one  to  either  extreme  one  is  equal  to  the  kinetic  energy  at  the  neutral 
configuration.  (The  absence  of  inertial  coupling  and  the  cancellation 
between  the  elastic  and  centrifugal  couplings  mean  that  there  is,  in 
effect,  no  coupling  when  the  blade  is  in  vacuum.  ) 

As  an  example,  consider  the  case  with  the  following  root 
conditions:  ^  - 

d  r 


~  .  dc 

L°  Ur 


h  -  C  d  j 

where  ^,-and  c^arc  given  torsional  elastic  constants.  For  this  case, 
( 1  03)  and  (1  1  0)  yield 


and 


A  =  A-  *  0 

tj  J  *• 


a 


* j 


ir  dr  r  =  r„ 


dct  dc, 


dr  d>' 


'  'r  ra 


and,  hence,  (112)  yields 
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Ij 

u.  h  r  -  c  C  Jr 
*  J  *  j 

o .  j  * ;  . 

Consequently, 

(106)  yields 

dh  d  h  ) 

L  (  *  J  .  1 

i  dCc  ^  ' 

\  dr  air  r  e  r0 

1  (Jr  dr  r  t', 

-  ^  o . 

X 

rv.. 

(113) 


(IN) 


and  also 


A 


Z  ” 

n  rii 


(115) 

For  the  two  extreme  cases  where  the  blade  is  cantilevered  at  r  -  rQ 
(namely,  both  tF  and  JcD  approach  infinity,  yet  both  /♦*  and  A7  are 
finite)  and  where  the  blade  is  fully  articulate  (namely,  -  0) 

the  relations 


C 


dh  \ 


dr  dr  r^r 
hold  for  j  ^  i  or  J 


(<*CC 


dc 


-  0 


~  i 


\  dr  dr /  ,  .  , 
o  o 

so  that  ( 1  1  4)  and  (115)  reduce  to 


-nr 


and 


=•  0 


j  y 


k.  -  n 


z  A 
cv  M 


L  L 


(116) 


(117) 


Results  in  agreement  with  (113)  and  (1  16)  are  obtained  tn  Reference  6, 
which  treats  a  blade  cantilevered  at  r  -  0. 

As  another  example,  consider  a  blade  with  flapping  and  lead- 
lag  hinges.  With  reference  to  Figure  33,  the  flapping  hinge  is 
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} , 


Figure  33.  RIGID  LINK  CONNECTING  FLAPPING  AND  LEAD-LAG  HINGES. 


located  at  r  =  />  and  the  lead-lag  hinge  at  r  =  rQ  .  The  link,  shown 
in  Figure  33,  connecting  the  two  hinges  is  assumed  to  be  rigid,  so 
that 


(h) 


r‘ro 


(118) 


and 


(‘Ktr.  'Wr.r.  = 


(119) 


are  three  of  the  root  conditions.  The  other  root  condition  can  be 
obtained  from  the  requirement  that  the  moment  about  the  flapping 
axis  must  vanish,  namely,  (M)  =0.  This  moment  may  be 

expressed  as 


where  is  the  moment  about  the  flapping  hinge  due  to  the  inertia  of 
the  link.  Letting  JH  be  the  mass  moment  of  inertia  of  the  link  about 
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the  (lapping  lunge  uiuJ  he  a  mass  element  in  the  link, 


U 

n 

uJVw' 

*1 

l< 

-  r 

<  ■  , 

or, 

a  s  s  uming 

the  link 

to  he 

that 

J  >'< 

>•),,»> 

-  0  , 

l  4 

-  ,  y  i  *■  T 

d>) 

-  il^i 

(JA ) 

L 

*(  dr. 

f  r, 

C ' 

\dr)rj 

or, 

using 

J, 

/(<’'- 

>»*  + 

“  |  UJ 

Z  T 

^  '  H 

whe 

re  mH 

is  til 

e  mass 

of  the 

h  J>  rt  j 

n(r  'v  /_  '  > ! 


>' H  e H 


)i  ±1 

'  V  dr  r  r 


(120) 


from  the  flapping  hinge  to  the  c  .  g.  of  the  link.  Therefore,  the  fourth 
root  condition  is 


[ A?  -  (rc  -  r,  )  j  ^  ^  -  {-  J %  -  ci ‘(j*  -  ?  ^  mH  r„  e„)\(jX  .  r<  _ 


ah 


(121) 


By  the  use  of  (1  18),  (110)  and  (12  1  ).  (100)  and  (110)  yield 
in  in  , 


a  nd 


<v 


=  +  CL*  {J*  -  *  mH  rF  e.„  )  j  J  -hl  aHj 


■ir  d> 


r  ^  x  r. 


and,  hence,  (M2)  yields 


^  in  Ih 

V  =  Af  4 

■J  V  *  l,  ir 

J  /  ‘  • 

>’  -  >\~ 

(122) 

ConsequentK  .  ( 1  OS)  yields 

V 

<■ J 

-  rr  t  -  o .  j  i  c 

-j 

(123) 

a  nd 

< 

L  i. 

-n2r  =  a/  m  ■ 

Li  L  i  L 

(124) 
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w  h  e  r  c 


A 

T 


•If' 


o ) 


{trdr  *-  rF  e„ 


rc 


(  dh-  dh4  ) 

J  ■  Ur  aLr  / F  . 


Finally,  some  remarks  on  the  teetering  rotor  (in  which  the 
two  blades  can  flap  on'y  as  a  continuous  beam  and  rf ■  -  0)  will  be 
made.  The  hub  mechanism  that  flaps  with  the  blades  i&  usually 
treated  as  a  rigid  mass  rigidly  attached  to  the  blades  at  r  -  0.  This 
rigid  mass  is  symmetric  with  respect  to  the  ^  -ax»s,  just  like  the  two 
blades.  Suppose  each  blade  has  a  lead-lag  hinge  at  r  =  r0  (the  results 
below  apply  to  the  case  w'ithout  lead-lag  hinges  upon  setting  rr  =  0). 
Assume  that  the  portion  of  the  blades  between  the  two  lead-lag  hinges 
is  rigid,  and  let  Z  be  the  total  moment  of  inertia  about  the  flapping 
axis  of  the  masses  of  the  hub  mechanism  and  the  part  of  the  blades 
between  the  lead-lag  hinges,  with  Z  being  contributed  by  the  j- 
direction  spread  only.  Using  the  superscripts  [  1  ]  and  [Z  ]  to  indicate 
quantities  belonging  to  the  two  different  blades,  the  joining  together  of 
the  two  blades  at  r  =  0  requires 


(dh  )  W  _fdh)in 
\dr)rT  ‘  ‘Ur/ 


Consequently,  the  two  blades  can  have  only  two  different  types  of  free 
oscillations: 


1  ype  I: 


A'1  =  -c-t'J  ; 


Type  II: 


I*-'-#)™ 


M)  [n 

r*  ru  i  “  r  7  r  -  ra 

[  Mh  is  given  by  (lZO)with  rF  =  O.  ] 


-  mi". 


hul  -  /T  ”  , 


(%) 


C?J 


r  - 


dh 

dr 


c 

VJ 

r  =  rn 


C/3, 


—  0. 
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The  results  turn  out  to  be  that,  for  both  i  and  j  belonging  to  either 
Type  I  or  Type  II,  (1  22)  through  (124)  are  satisfied  by  either  blade. 
Note  that  the  last  condition  of  1>  pc  II  automatically  reduces  (122), 
(123)  and  (124)  to  (113),  (116)  and  (l  IT),  respectively.  Il  one  of  the  i 
and  j  belongs  to  Type  I  and  the  other  to  Type  i.,  then  define 


K 

-  < r" 

+ 

K 

LJ 

‘■j 

‘■J 

*  10 

*  U1 

T 

=  T 

+- 

T 

<■1 

Lj 

A 

•  % 

* 

M 
‘ J 

for  the  two  blades  together  and  find 


so  tha f  (122)  and 


% 


A 

z  M 


.  (dhL  dhJ  ) 

l  dr  d  -  ) r*r 


=  0 


CJ  \  dr 

(123)  are  satisfied  by  the  two  blades  together. 
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During  a  previous  program,  Cornell  Aeronautical  Laboratory,  Inc.,  devel¬ 
oped  a  method  ci  computing  rotor-blade  loads  and  motions  of  a  single-rotor  helicopter 
in  steady  forward  flight  by  assuming  that  the  blade  pitch-control  settings  (collective, 
longitudinal  cyclic,  and  lateral  cyclic)  and  the  rotor-shaft  tilt  angle  are  known  and 
that  the  blade  motions  are  restricted  to  flapping  and  flapwise  bending.  >The  present 
effort  was  undertaken  to  extend  the  previously  developed  method  by  (1)  including  the 
blade  inplane  and  torsional  motions,  and  (2)  treating  the  four  trim  constants  (namely, 
the  blade  pitch-control  settings  and  the  rotor-shaft  tilt  angle)  as  unknowns.  The  trim 
constants  are  different  for  different  flight  conditions  and  are  determined  through  the 
use  of  four  appropriate,  average  equilibrium  conditions  of  the  helicopter.  These 
equilibrium  conditions  are  called  the  trim  equations  and  are  derived  in  this  report, 
taking  into  account  the  inertial  forces  of  the  blades  due  to  elastic  deformations. 

Lagrange's  equations  for  the  blade  motions  are  given.  The  generalized  coordinates 
employed  in  these  equations  to  represent  blade  bending  are  those  which  have  oeflection 
components  in  two  mutually  perpendicular  directions.  Orthogonality  relations  of 
vibration  modes  of  twisted,  rotating  blades  are  derived  and  are  used  in  simplifying  the 
equations  of  blade  motion.  A  successive  approximation  procedure  was  developed 
which,  upon  incorporation  with  the  previously  developed  iterative  procedure,  yields 
the  aerodynamic  loads,  the  blade  responses,  and  the  required  trim  constants. 

Computed  results  are  obtained  for  the  UH-1A  rotor  at  advance  ratios  of  0.  26  and  0.  08 
and  for  the  H-34  rotor  at  advance  ratios  of  0.29  and  0.  18.  Comparisons  of  these 
results  with  available  measured  results  are  presented. 
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